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Abstract 



We study both the static and dynamic properties of gapped, one-dimensional, 
Heisenberg, anti-ferromagnetic, spin chains at finite temperature through an 
analysis of the 0(3) non-linear sigma model. Exploiting the integrability of 
this theory, we are able to compute an exact low temperature expansion of 
the finite temperature correlators. We do so using a truncated 'form-factor' 
expansion and so provide evidence that this technique can be successfully 
extended to finite temperature. As a direct test, we compute the static zero- 
field susceptibility and obtain an exact match to the susceptibility derived 
from the low temperature expansion of the exact free energy. We also study 
transport properties, computing both the spin conductance and the NMR- 
relaxation rate, 1/Ti. We find these quantities to show ballistic behaviour. In 
particular, the computed spin conductance exhibits a non-zero Drude weight 
at finite temperature and zero applied field. The physics thus described differs 
from the spin diffusion reported by Takigawa et al. Q from experiments on 
the Haldane gap material, AgVP2SQ. 
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I. INTRODUCTION 



The realization that one-dimensional, integer spin, antiferromagnets possess an energy 
gap has made these systems the object of intense study. The model perhaps most com- 
monly used to explore their properties is the field theoretic 0(3) non-linear sigma model 
(NLSM) Although the model has the virtue of being integrable P,p!0|, its properties 

are nonetheless only partially understood. It is possible to access static, thermodynamic 
quantities while dynamic properties, in particular, transport properties, are in general, un- 
available. These latter quantities depend upon knowledge of correlation functions which are 
generically not exactly computable in integrable models. There are, of course, perturba- 
tive techniques by which correlators in the 0(3) NLSM may be analyzed. But in strongly 
coupled models, of which the 0(3) NLSM is one, perturbative techniques present a host of 
difficulties and so can miss qualitative (never mind quantitative) features in the physics. 

The inability to completely understand correlation functions in the fully quantum 0(3) 
NLSM has been at the root of a recent controversy in the literature. Takigawa et al. |1| 
demonstrated through measurements of the NMR relaxation rate, 1/Ti, of the Haldane 
gap compound, AgVP2SQ, that at long wavelengths, the spin-spin correlation functions are 
diffusive in nature. In an elegant series of papers, Sachdev and Damle [|^,|| developed 
a semi-classical treatment to attack the problem and subsequently were able to describe 
this diffusive behaviour. Nonetheless their computation was semi-classical leaving open the 
possibility that a fully quantum treatment of the 0(3) NLSM would lead to different physics. 
This possibility was hinted at in the work of Fujimoto [|ll|]. There the spin conductance 
was computed using exact thermodynamic considerations. Upon subsequent work |T^, it 
became clear that the two treatments produced qualitatively different results. In particular, 
the Drude weight of the spin conductance, D, of the 0(3) NLSM was found to be non- 
vanishing in the zero field limit [|1^], whereas the corresponding semi-classical treatment sees 
D{H = 0) = 0. This qualitative difference opens up the possibility that diffusive physics 
is not present in the 0(3) NLSM itself but requires some additional mechanism. Such 
mechanisms might include a spin-phonon coupling (as suggested by [|lT|), spin anisotropy, 
inter-chain coupling (the spin-chains in AgVP2SQ are only quasi 1-D - there do exist weak 
couplings in between chains although the weakness of these couplings seems to preclude this 
possibility), or perhaps small generic integrable-breaking perturbations [^. 

In this paper we attempt to address this problem by demonstrating a technique to com- 
pute exactly a low temperature expansion of correlators in the 0(3) NLSM. This expansion 
is based upon a 'form-factor' expansion. Form-factor expansions have a long history in the 



computation of correlators |13-Tq]. However these expansions have been used almost exclu- 



sively at zero temperature. When they have been used at finite temperature, they have been 



used either in the computation of expectation values lacking dynamical properties p9H22 



or in the development of distinct non-perturbative representations (i.e. Fredholm determi- 
nants) of correlators p3|j2^ where all the terms in the expansion were kept. In this article 
we show that truncated form-factor expansions can be used to sensibly describe correlation 



functions at finite temperature. This is distinct from the programme proposed in p5|j26 
where form factor expansions were employed but the form factors themselves were recom- 
puted to take direct account of thermal fiuctuations. Here we employ the same form factors 
used in zero temperature computations. 
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A form-factor expansion of a correlation function is predicated upon some generic prop- 
erties of integrable models. Most importantly, the exact eigenfunctions of the model's fully 
interacting Hamiltonian are known. With this knowledge comes a well-defined notion of 
'particles' or elementary excitations in the system. The scattering of these particles is com- 
pletely described by two-body S-matrices. In particular, particle non- conserving processes 
are disallowed. Ultimately this feature is a consequence of a series of non-trivial conservation 
laws possessed by the integrable model. In some sense, an integrable model is a superior 
version of a Fermi liquid: a particle's lifetime is infinite regardless of distance from the Fermi 
surface. 

In order to understand these features of the 0(3) NLSM, we begin by providing an 
overview of the model. The 0(3) NLSM is described by the action, 

S= — [ dxdtid^'nd.n), (1.1) 
2g J 

where n = (ra^., n^, n^) is a bosonic vector field constrained to live on the unit sphere. This 
action is arrived at from the Hamiltonian of the spin chain, 

H = JJ2S,-S,+,. (1.2) 

i 

In the continuum, large s, limit, the spin operator. Si, is related to the field, n, via 

= i-iysn, + Mi, 

that is, n(a;, t) is the sub-lattice or Neel order parameter. M on the other hand describes 
the uniform (i.e. wavevector A; ~ 0) magnetization. M is related to n via 

M = -n X dtn, 
9 

and so is given in terms of the momentum conjugate to n. 

The low energy excitations in the 0(3) NLSM take the form of a triplet of bosons. The 
bosons have a relativistic dispersion relation given by 

Here A is the energy gap or mass of the bosons related to the bare coupling, g, via A ~ 
jg-27r/g rjj_^g dispersion relation of all three bosons is identical as the model has a global 
SU(2) symmetry. The exact eigenfunctions of the 0(3) NLSM Hamiltonian are then multi- 
particle states made up of mixtures of the three bosons. Scattering between the bosons is 
described by the S-matrix 

2me 



n m = - . 

' (0 + m)(0-^27r)^ 
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- {e + ^r.){e-^2^)■ ^^'^^ 

Here 9 parameterizes a particle's energy/momentum via E = Acosh(0), P = Asinh(0). 
The primary advantage of this parameterization is the implementation of Lorentz boosts. 
Under such a boost, 6 6 + a. As such Lorentz invariant quantities are invariably functions 
of differences of rapidities. We stress that this relativistic invariance is a natural feature of 
the low energy structure of the spin chain. (However we do point out for spin 1 chains, 
A ~ AJ. As J serves as the cutoff for the theory, the low energy sector of the theory is not 
unambiguously defined.) 

With the excitation spectrum of the 0(3) NLSM in hand, we return to the form-factor 
expansion. A finite temperature expansion of correlators is given in terms of a trace over 
the Boltzmann density matrix: 

G^{x,t) = lTr(e-^^O(x,t)C(0,0)) 

^ e-^^-- (n, Sn\0{x, t)0(O, 0) |n, g„) 

Here the state, \n,Sn), denotes a set of n-particles carrying spin quantum numbers 
Inserting a resolution of the identity between the two field then leads us to a double sum, 

E"»« e-^^'"{n,Sn\O{x,t)\m,Sm){m,Sm\O{0,0)\n,Sn) 
G^(x,t) = ^^ — , ^ . (1.5) 

We thus have reduced the evaluation of the correlator to the evaluation of a series of matrix 
elements (known as 'form factors'). In an integrable model like the 0(3) NLSM, these matrix 
elements are in principle exactly computable. However as the number of excitations involved 
increases, the functional forms of the matrix elements become increasingly unwieldy. This, 
together with the difficultly in evaluating the sums, J2{n,s„),{m,sm)^ ensure in all but a few 
special cases the correlators do not admit a closed form expression. 

To surmount this we adapt an idea from zero temperature form-factor expansions. Rather 
than look at the correlator in real space and time, we examine the (more relevant) related 
spectral function, G'~^{k, uj). In computing G'^(fc, tu), only terms in the form factor sum with 
a given energy, tu, and momentum, fc, contribute to the sum, 

nsn 
msm 

xe-^^- (n, s„|O(0, 0)|m, s^)(m, s™|O(0, 0)|n, s,), (1.6) 

as enforced by the presence of the two delta functions. For any cj, fc, this dramatically 
reduces the number of matrix elements one must compute^. This reduction nonetheless 



^Here is simply the Fourier transform of G'^(x, t), but similar considerations also apply to the 
corresponding retarded correlator. 



4 



leaves a difficult computation. However we can exploit the gapped nature of the spin chain 
to make the problem more tractable. Because the theory is gapped (with gap, A), the 
correlator admits a low temperature expansion of the form, 

G''ik,uj) = J2a^{k,u)e-^^. (1.7) 

n 

For the particular correlators of concern in this paper and for the range of u and k in 
which we are interested, each a„ is determined by a single matrix element. Because we can 
compute these matrix elements, we obtain an exact low temperature expansion. 

Our ability to compute such an expansion bears upon another controversy in the liter- 



ature. LeClair and Mussardo argued that it was possible to use the same form-factors 
we employ here to compute finite temperature correlators. However rather than directly 
evaluate individual terms in the sum (|1.5|) , they first conjectured an ansatz involving a re- 
summation of terms in the sum. This is described in more detail in Section 3. This procedure 
was criticized in [^]. There it was argued that while this worked for the computation of 
one-point functions, it was problematic for two-point functions. Rather it was argued it was 
better in general to attack such problems through the use of form factors computed against 



a thermalized vacuum |p5H27[|. However the counterexample cited in |^0|) ^ computation 



involving interacting quantum Hall edge states, involved a gapless theory, and so is in a dif- 
ferent class than the model considered in this paper. (Without a gap, the low temperature 
expansion we consider above ceases to make sense.) This work here shows that it is possible, 
at least in certain cases, to make sense of the form-factor expansion of two point functions 
at finite temperature. But while we can make sense of this expansion, we cannot compare 



our computations directly to the ansatz posited in [|T9|. Their ansatz as is applies only to 
diagonal theories where scattering does not permute internal quantum numbers, contrary to 
the case here. 

The outline of the paper is as follows. In Section 2 we summarize the results of the form 
factor computations for three quantities: the magnetic susceptibility, and two transport 
properties, the spin conductance and the NMR relaxation rate, 1/Ti. The details of these 
computations are found in later sections or in appendices if the reader is so interested. The 
first quantity, the susceptibility, is compared to the susceptibility as derived from a low 
temperature expansion of the exact free energy. We see that they match verifying our claim 
that the form factor expansion can yield an exact low temperature expansion. 

We compare our transport calculations to the semi-classical computations in [|^|[. The 
essence of this method lies in treating the spin-chain as a Maxwell-Boltzmann gas of spins 
which interact with one another through the low energy limit of the scattering of the 0(3) 
NLSM, 

S:tiO = 0) = -SadScb. (1.8) 

While static properties computed in the two treatments agree (for the susceptibility, we find 
that up to temperatures on the order of the gap, T ~ A, the two computations agree), we see 
differences in transport properties. For the spin conductance we find, in contradistinction 
to the semi-classical computation, that the Drude weight of the spin conductance is finite 
in the limit of zero external field. Our results for the NMR relaxation rate, 1/Ti, indicate 
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a similar discrepancy. We, like 0, find that 1/Ti is characterized by ballistic logarithms. 
These logarithms are relatively robust: they continue to appear at higher orders in the low 
temperature expansion. We do not, however, see diffusive behaviour in the relaxation rate, 
i.e. l/Ti ~ nor does our low temperature expansion match the low temperature 

expansion of the semi-classical computation of the correlator. 

We consider two possibilities in explaining these discrepancies. We argue that the struc- 
ture of the conserved quantities or charges differs between the 0(3) NLSM and its semi- 
classical variant and that these differences lead to ballistic behaviour on the one hand and 
diffusive behaviour on the other. The other explanation we forward to explain this discrep- 
ancy lies in the supra-universality of the low energy S-matrix ( |1.8| ). The low energy limit of 
this S-matrix is shared by generic integer spin chains. Indeed it is shown in ^ that a two-leg 
spin-1/2 ladder, expected to share the low energy behaviour of a spin-1 chain, has this exact 
low-energy S-matrix. However rather than the supra-universality being a virtue, it may be 
that it under-specifies the physics. In this way the semi-classical treatment, valid in and 
of itself (particularly in light of its ability to reproduce experimental data), may capture 
different physics than that of the 0(3) NLSM. In Section 2 we consider this further in the 
light of the sine-Gordon model where a similar phenomena may be argued to occur. 

In the first part of Section 3 we explain in some detail how the form-factor expansion is 
to be understood. In particular we consider the various technical details of the expansion, 
including how to regulate the infinities that appear generically in the form factors of the 
double expansion. In second part of Section 3 we review the specific form factors of the 0(3) 
NLSM. And finally in Section 4, we review the low temperature expansion of the exact free 
energy, necessary for comparison with the form-factor computation of the susceptibility. 
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II. SUMMARY AND DISCUSSION OF RESULTS 



A. Zero Field, Finite Temperature Susceptibility 

In this subsection we present results for the magnetic susceptibility arising from several 
methods of computation: a form factor evaluation of the magnetization-magnetization cor- 
relator in the context of a Kubo formula, an exact computation of the system's free energy, 
and finally, treating the excitations of the 0(3) sigma model as non- interacting particles 
obeying both a Fermi-Dirac distribution and a Maxwell-Boltzmann distribution in the spirit 
of the semi-classical approximation of Sachdev and Damle [^H- We thus will be able to 
determine the temperature regime over which our truncation of the form factor expansion 
applies as well as comparing with other computational techniques. 



1. Kubo Formula and Form Factors 



The susceptibility, at = can be computed from the magnetization-magnetization 
operator using a Kubo formula: 



XiH = 0) = C{uj = 0,k = 0) 

rl3 



C{uj = 0,k = 0) 



dx 



{nM^ix,T)M^{0,0))) 



(2.1) 



To evaluate this correlator we employ an expansion in terms of the exact eigenfunctions of 
the theory, i.e. a form factor expansion. In particular we write 



1 

Z 



(M3(x,r)M^(0,0)) = ^Tr(e-''5^C(x,t)O(0,0)) 

e-f'^'-{n,Sn\O{x,t)O{0,0)\n,Sn 



(2.2) 



Here \n, Sn) is a state of n excitations with spins described by Sn = {si, ■ ■ ■ , Sn}- In writing 
the above we have suppressed sums over the energy and momenta of the excitations. A term 
in the thermal trace with n excitations is weighted by a factor of e~"^^. At low temperatures 
it is thus a good approximation to truncate this trace. For this computation we keep only 
terms with one and two excitations, i.e. n = 1, 2. To evaluate the matrix elements appearing 
in ( |2.2D we insert a resolution of the identity in between the two fields. As we only consider 
matrix elements involving one and two excitations from the thermal trace, we thus have 

{si\M^{x,T)M^{0,0)\si) = J2 {si\M^{x,T)\mSm){mS^\M^{0,0)\si) 

mSm 

= Y: {silMlix, t)\s[){s[\MI{0, 0)\s,) + ■■■■ 
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{siS2\Ml{x,r)Ml{0,0)\s2Si) = Yl {siS2\Ml{x,r)\mSm){mSm\Ml{0,0)\s2Si) 

mSm 

= J2 {s,S2\M^{x,t)\s[s'^){s'^s[\M^{0,0)\s2S,) + ■■■. (2.3) 

In the above we have truncated the sum arising from the resolution of the identity. With 
the first matrix element of the thermal trace, we only keep terms from the resolution of 
identity with one excitation. We are interested in the behaviour of the susceptibility at 
a; = and this term provides the only contribution. Similarly, the only term arising from 
the second matrix element of the thermal trace contributing to the DC susceptibility comes 
from keeping the term from the resolution of the identity involving two excitations. Further 
details surrounding the methodology of this expansion and the explicit exact evaluation 
of the matrix elements are found in Section 3 and Appendix A. With such details we can 
evaluate C{u; — 0,k — 0) with the result 

C{uj = 0,k = 0) = Ci{uj = 0, = 0) + C2{uJ = Q,k = 0), (2.4) 

where Ci and C2 are given by 

C^(cu = 0, A; = 0) = ^^Ki{(3A); 

TT 

C2{uj = 0,k = 0) = -^Xi(2/?A) 

TT 



^^2 + 5712^12 + 471^' 



/ d^id^2e-^^('^°^^(^^)+^°^'^(^^»cosh(^0 

TT J 

■^iri(2/3A) + ^Xo(/3A)Xi(/3A) + 0{^e-'P% (2.5) 



where 6*12 = Oi — O2 and Kn are standard modified Bessel functions. The first term in C2 
is a 'disconnected' contribution related to Ci. The second term is a connected contribution 
and as such is genuinely distinct from Ci. We now consider such disconnected contributions 
further. 



2. Resummed Form Factors 



In computing the susceptibility, we are able to go beyond the approximation introduced 
in truncating the form factor sum arising from the thermal trace. It is possible to include 
'disconnected' terms arising from higher particle contributions. Such disconnected terms 
appear when higher particle matrix elements are evaluated. For example when we evaluate 
the four excitation matrix element (s2s'i|Mq (0, 0)|s2Si), we obtain a term of the form 

(44|m3(o, o)|s2Si) = • • • + 5,.,,(s;|m3(o, o)|s2) + • • • . 

This term is 'disconnected' in that it is directly related to a matrix element involving a lesser 
number (two) of excitations. It arises from the annihilation of s'2 with Si. Such a term is 
responsible, as we just indicated, for the first term of C2 above. 
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What is remarkable is that we are able to sum up over all possible disconnected pieces 
arising from arbitrarily high particle form factors which are proportional to the connected 
lower particle matrix elements already computed. This resummation amounts to the eval- 
uation of a geometric series. For example, including all disconnected terms involving the 
matrix element going into the evaluation of Ci modifies it as follows 



C, = (^ f ci^e-^^^°^^(^) cosh(e) ^ ^ / rf^e-^^^^i^W cosh(^) V S^e""^^ 



cosh(6») 

n=0 



_ ^ r e-P^^oMe) cosh(g) 

- ~ I l + 3e-/3Aco.h(e) 

= ^ir^(^A) - ^1^K,{2(3A) + 0(e-3/3^). (2.6) 

We see resumming the disconnected pieces thus reproduces both Ci and the first term in C2 
plus additional terms higher order in e~^^. The appearance of the factors e-^^=°^^(^) in the 
geometric series is natural and arises from the Boltzmann weighting of the higher particle 
terms. The combinatorial factor of 3 reflects the three bosons in the system. 

In collecting all the disconnected pieces related to the connected term in C2, we find 
something similar 

connected C2 + disconnected terms 

[ g-/3A(cosh{0i)+cosh(92)) ll7r^ + 26'f. 



TT 



iO^iO^^^^p^K^M) + co=h(«.))-— (2.7) 



One expects in general that the inclusion of disconnected terms from arbitrarily high particle 
number will improve the accuracy of the calculation. In the case of the susceptibility, the 
agreement between the form factor computation and the exact numerical analysis actually 
becomes slightly worse. However this should not be taken as indicative of the resummation 
in general. We will comment on this further at the end of this section. 



3. Gases of Free Particles 

For the purposes of comparison, we compute the susceptibility of both a free electron 
gas (or equivalently, a system of hard-core bosons) as well as a Maxwell-Boltzmann gas. At 
sufficiently low temperatures both of these quantities should be close to the exact value of 
X for the 0(3) sigma model. How the susceptibility of the free electron gas deviates from 
the exact value of x gives us an understanding of the temperature at which interactions 
become important. And how the susceptibility of the Maxwellian gas deviates from the 
exact answer marks the temperature at which the semi-classical approximation found in 
Damle and Sachdev must begin to breakdown. 

These two susceptibilities are given by 



_f3A /• cosh(^)e-^^'=°^^(^) 

Xfree el. - — j "-^ (^l g-/3Acosh(e))2 
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IT 



TT A 



XMB = f-^e-^^. (2.8) 



vr 



We see that at low temperatures (/3A <C 1) both of these expressions coincide with the low 
temperature limit of the form factor computation of x- Iii particular, the terms of 0{e~^'^) 
are identical. 



4- Thermodynamic Bethe Ansatz 

It is possible in the case of the 0(3) sigma model to arrive at exact expressions (in the 
form of coupled integral equations) for the zero-field susceptibility [0,0 . These equations, 
in their most compact form, appear as 



^^H = 0) = -— J d9cosH9) 



H=0 



e(^) = Acosh(0) -T J ci^'log(l + e^'^'^^'^)s{e - 6'); 
eni9) =T j de's{e - e'){ log(l + e"'-'^'"^) + log(l + e^^"+^(^')) + d^n log(l + e^^'"^] 

H = hm ^ (2.9) 

n—too fi 

We will show results from the exact numerical evaluation of these equations in the next 
section. However these equations admit a closed form low temperature expansion. The 
details of this expansion may be found in Section 4. Here we just give the final results 

TT TT 

, r -flAl'r.oshl'fili Wcosh('(9-,11 i /n \ 1171^ + 26'?; 



TT 



Remarkably, we see this expansion agrees exactly with the corresponding expression derived 
with the aid of form factors. Thus the form factor expansion at finite temperature meets an 
important test. 



5. Comparison of Methodologies 

In this section we compare the various methods of computing the susceptibility of the 
0(3) sigma model. In Figure 1 are plotted the susceptibilities computed via an exact numer- 
ical analysis of the TEA equations, a low temperature expansion of the same equations, and 
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a computation based upon the two and four particle form factors. We see that as indicated 
previously that the form factor computation and the low temperature expansion match ex- 
actly. Moreover these two computations track the exact susceptibility over a considerable 
range of temperatures despite the fact these computations are truncated low temperature 
expansions. 

In Figure 2 we compare both the exact TBA susceptibility and the form factor computa- 
tion of X with the susceptibility of a classical Maxwellian gas. We see the results track one 
another for temperatures, T < A. For temperatures beyond A, however, the Maxwellian 
susceptibility differs markedly. This is then roughly the temperature at which the semi- 
classical approximation found in [|^|[ should be expected to break down. 

In Figure 3 are plotted the exact results for the susceptibility together with the sus- 
ceptibility from the resummed form factors. We see that the resummed susceptibility is 
somewhat higher that the exact numerics at T ~ 5A and disagrees at roughly the 10% level 
whereas the susceptibility computed using the unresummed form factors sees better agree- 
ment at these same temperatures. At lower temperatures (T ~ 2 — 3A) the disagreement 
between the exact numerics and the two form factor computations is roughly the same. In 
general then, the resummation does not improve the accuracy of the computation of the 
susceptibility. 



B. Spin Conductance 

In this section we compute the spin conductivity, cr^. The spin conductivity gives the 
response of the spin chain to a spatially varying magnetic field. It is defined via 

3i{x,t) = <JsVH, (2.11) 

and so can be expressed in terms of a Kubo formula, 

Rea,(fc,cu) = -i| Me^'=^+^-*ImOo(x,t)ji(0,0)),,,,,,,,. (2.12) 

In the notation used in this paper the spin current ji is synonymous with Mi, the Lorentz 
current counterpart of the uniform magnetization, Mq = jo- We will focus primarily on 
computing the Drude weight, D, of Reds, i.e. computing the term in as{k,uj) of the form 

as{k = Q,uj) = D5{uj). (2.13) 

However we are able to compute cxs for general k^uo. We find that for uj -C 2A, A; = 0, 
the spin conductivity is described solely by the Drude weight. In particular, we find no 
indication of a regular contribution to (Jsik = 0,ct;). 

To evaluate ds, we employ the identical form factor expansion to that used in comput- 
ing the susceptibility. And like the susceptibility, our result is an exact low temperature 
expansion of -D, 
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Here we will compute Di and D2 exactly. As the details of the computation are nearly 
identical to that of the susceptibility, we merely write down the results: 

D(H = 0) = /3A / rf^e-^^^°^^(^)^^^^i^(l - 3e-^^^°^^(^)) 
J cosh(6') 

+2dA ( de.de.,c-'^^^^°'^^'''^^'°'^^^''^^ smh^(g2) llTT^ + 261^ + Ofp-^^^V 



-13 A, 



-e-''\ 4^ilV^ - X + + Oie-'^''). (2.14) 



" TT V A 'A 



This expression involves only the two and four particle form factors. If we also include 
all higher order disconnected terms related to those above we find instead (akin to the 
susceptibility) , 



D{H = 0) = /3A / dee-^^'^^'^^'^^^^^^- 
J coshifc') 1 



ih{e) i + 3e/5Acosh(e) 

oiv.T.S/'/a ^ 1 1 ^2 I 0/32 

+2/3A / d9id92e 



cosh (^2) Of 2 + 57i^9j2 + 47r4 
11 1 
^21 + 3e/3^c°'^Mfi) 1 + 3e/3Acosh{02)''' (2.15) 

We plot these two results in Figure 4 as a function of T/A. Akin to the susceptibility, the 
result does not differ greatly if the resummed disconnected terms are included. 



We observe that D(H = 0) 7^ 0. This is in accordance with |]TT| where D is computed 
using an argument involving the finite size scaling of the thermodynamic Bethe ansatz 
equations. (We do note that the computation of D at = in |jll| appears only as a note 
added in proof and so is decidedly sketchy. However the equations governing D developed 
in IjTlJ are manifestly positive with the consequence D cannot vanish.) But our results do 



differ from the semi-classical computation of [O where is was found that D vanishes at 



H = 0. And again we find no additional regular contributions to as{k = 0,a; = 0) - only 
the Drude term is present. This is true not just to the order of the form factor expansion 
we work but at least to one higher order. Moreover we are willing to conjecture that is true 
to all orders. 

We have only given the spin conductivity at if = 0. However it is extremely straightfor- 
ward to generalize the form factor computation to finite H. As H couples to the total spin, 
a conserved quantity, the form factors, f^{x,t), of an operator, 0{x,t), carrying spin s, are 
altered via the rule 

(In the case of the spin conductance, the spin currents, = M^, carry no spin and so are 
not altered at all.) The only remaining change induced by a finite field is to the Boltzmann 
factor appearing in the thermal trace. If an excitation with rapidity, 6, carries spin s, its 
Boltzmann factor becomes 

-/3(Acosh{e)-sH) 
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For example we find D as a function of H (to 0{e ''^)) to be 

,sinli2(^) 
cosli(^^) 



D(H) = f3 A coshipH) J d^e-^^™^^(^) 



(2.16) 



Again this in agreement with [|rT|]. Indeed |[TTI| computes D{H) at large H/T (but H A) 
to be 



471 7 cosh (61) 



(2.17) 



Up to a factor of 27r, this expression is in exact agreement with 2.16 . In this particular case 
our derivation of D{H) agrees with the semi-classical computation [|12| (provided T <^ H <^ 
A). The symmetries in the semi-classical model that lead D{H = 0) to vanish are broken 
for finite H. 



C. NMR Correlators 

In this section we compute the NMR relaxation rate, 1/Ti. We are interested in comput- 
ing this rate in order to compare it to the experimental data found in on the relaxation 
rate of the quasi one-dimensional spin chain, AgVP2SQ. For temperatures in excess of lOOi^' 
(the gap. A, in this compound is on the order of 320K), the experimental data IQ shows 
the relaxation rate to have an inverse dependence upon y/H: 

1/Tioc 



H 

This dependence is nicely reproduced by the semi-classical methodology in Moreover 
the semi-classical computation reproduces the activated behaviour of 1/Ti in this same 
temperature regime: 

l/Ti oc 6-3/3^/^ 

We are interested in determining whether a calculation in the fully quantum 0(3) NLSM 
can reproduce these results. To this end we compute 1/Ti using a form factor expansion. 
Sagi and Affleck have already done such a computation to lowest order in e~^^. But 
they do not find the above behaviour. Rather they see 

l/Ti oc hgiH); l/Ti oc e"^^. 

We continue this computation one further step, computing to (9(e~^^^). Given the be- 
haviour, 1/Ti ~ appears only as T is increased beyond lOOK (i.e. T/A ~ 1/3), it is 
not unreasonable to suppose higher order terms in a low temperature expansion of 1/Ti are 
needed to see this singularity. 

To proceed with the computation of 1/Ti, we review its constituent elements. 1/Ti can 
be expressed in terms of the spin-spin correlation function 0: 
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^lTi= E / 7^^a/3(A;)A,-,(-A:)(Mo^Mo^)(A;,o;^), (2.18) 

c=l,2 -I 
13=1,2,3 

where uj^ = JnH is the nuclear Lamour frequency with the nuclear gyromagnetic ratio 
and the A^js are the hyperfine coupling constants. In the above we assume H is aligned in 
the 3-direction. The above integral is dominated by values of k near Moreover in the 
relevant experiment, the hyperfine couplings are such that only (MqMq) contributes. Hence 

1/Tioc {M^M^){x = 0,ujn (2.19) 

We now proceed to compute (MqMq). 

To compute (MqMq), we again employ a form factor expansion. Akin to the compu- 
tation of the susceptibility and the spin conductance, this computation amounts to a low 
temperature expansion of (MqMq), 

where we are able to compute ai and We place the details of this computations in 
Appendix B, here merely quoting results: 

{M^M^Kx = 0,u; = 0) = (— e-^^log(f ) -7) - — e'^^^ log(f ) - 7) 

\ 71 ri vr rl 

+Ae-^^^(log(^) - 7)^(247r + il))(l + 0(iJ/T) + 0(r/A)), (2.20) 

where 7 = .577. . . is Euler's constant. We are interested in the regime if ^ T ^ A (the 
regime where it is expected spin diffusion produces singular behaviour in l/Ti). The terms 
that we have dropped do not affect this behaviour. In principle there is no difficulty in 
writing down the exact expression (to (!?(e~^^^)); it is merely unwieldy. This expression for 
1/Ti is plotted in Figure 5 for a variety of values of the ratio T/A. 

We see that we do not obtain the same behaviour as found in Going to the next 

order in 0{e~'^^^) produces a behaviour in 1/Ti as if ^ identical to the lower order 
computation of (9(e~^^): we again find a logarithmic behaviour consistent with ballistic 
transport. An alternative comparison we might make to the results of is to perform a 
low temperature expansion (in 0{e~^^)) of the semi-classical computation of (MqMq) (a; = 
0,0; = 0). Doing so by treating Te~^^/H as a small parameter, we find 

{M'MK^ = 0, c = 0) cx Ae-^^(log(^) - 7 + - + (^(e"'"^))- (2-21) 

We see that the low temperature expansion of the semi-classical result agrees to leading 
order with our computation but afterward differs. (We have already seen that this occurs 
with the computation of the susceptibility.) It possesses no term of (9(e~^^^). The next term 
rather appears at (9(e~^'^'^) and possesses a 1/ii^ divergence. That the small H behaviour 
should be 1/y/H does suggest the importance of summing up terms. But the lack of a term 
of (9(e~^'^^) in the semi-classical result nonetheless hints that the two results are genuinely 
different. 
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D. Discussion 



We have demonstrated that it is possible to compute exact low temperature expansions 
of correlators using form factors. Moreover we have done so in a non-trivial theory where 
particle scattering sees the exchange of quantum numbers. An important question to answer 
concerns the breadth of the applicability of our techniques. Our ability to carry out these 
computations was partially predicated upon the particular correlators we studied. For exam- 
ple, the fact that only a single matrix element contributes at 0{e~^^) and 0{e~'^^'^) in the 
computation of the susceptibility is related to the magnetization operator in the 0(3) NLSM 
model being a Lorentz current density. Because of these particular details, we thus expect 
that exact low temperature expansions of correlators will not be available in all theories. 

The computation of correlators is done in the context of a grand canonical partition 
function. Specifically, we do not work at fixed particle number but include matrix elements 
involving an arbitrary number of particles or excitations (see for example). This differs 
from the treatment found in |2^. There correlators are computed in a canonical ensemble 
using form factors at some fixed particle number, N. A thermodynamic limit is then taken, 
N,L —y oo holding N/L (i.e. the particle density) fixed. On a technical level these methods 
may seem ostensibly different. In particular in this paper we end up computing correlators 



using form factors involving a small finite number of particles whereas computes cor- 
relators using form factors involving a diverging number of particles. It might appear then 
that we are somehow missing information that arises in working at a finite particle density. 
This would seem crucial in computing transport properties where a finite particle density is 
necessarily determinant. 

However this difference is only apparent. The N particle form factors used by include 



disconnected terms. These disconnected terms are equivalent to form factors involving small 
numbers of particles. The (large) N-particle form factors then contain the same information 
we use in our representation of the correlators. Moreover we can make this identification 
precise. Our use of form factors in the grand canonical ensemble involving some few number 
of particles, n, is predicated upon the small parameter, e""^^. But the disconnected terms 
of an N-particle form factor involving n particles (with n < N) are similarly weighted by 
the same small parameter, e"^^". More generally, the presence of a gap. A, thus means we 
can in principle create an explicit map between the two approaches. 

The semi-classical method found in is similar to the approach taken in in 



that it uses a canonical ensemble. It is an interesting question whether a grand canonical 
ensemble approach can be developed in this same semi-classical approach. The answer is 
not obvious. Our method works (at least at the technical level) because we can readily 
identify disconnected terms. It is not clear whether a similar identification can be made 
semi-classically. 

We do want to emphasize a caveat to our methodology as discussed in some detail 
in Section 3. It is unclear whether it is possible to compute quantities that show non- 
analyticities as T ^ 0. For example it is not obvious how to compute the thermal broadening 
present in the single particle spectral function. At T = it takes the form 



{nn){uj,k) ~5(cj- VA;2 + A2), (2.22) 
but is expected to broaden into a Gaussian-like peak at finite T. To see this in a form factor 
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expansion would likely require a resummation of terms. However it may well be feasible to 
deduce the necessary resummation from the lower order terms in the form factor expansion. 

We have also discussed using a resummation of higher order 'disconnected' terms to 
improve the form-factor computation. For the quantities considered, it turned out the 
resummation did not provide a real improvement to the original computation. Nevertheless 
we would guess that in general, the resummed form factors will provide a more reliable 
answer as the temperature is increased. It is an artefact of the above cases that they do 
not do so here. For example, we see that at extremely high temperatures, the susceptibility 
as computed by either of the form-factors methods saturates to a constant. As such, errors 
in either method are cutoff - as these expressions at T = oo do not differ greatly from 
their low T values, any potential error is bounded. If instead we computed the finite field 
magnetization where we would expect a linear T dependence, the differences between the 
two form-factor computations would be comparatively magnified. 

To come to some sort of judgement between the form-factor and the semi-classical ap- 
proaches, an understanding is needed of the differences between our computations of the 
spin conductance and the NMR relaxation rate. In the case of the first quantity, it is likely 
this difference is real and not an artefact of our methodology. The data that goes into the 
spin conductance is identical to that needed to compute the susceptibility and we know that 
we can match the low temperature expansion of the susceptibility with a similar expansion 
coming from the exact free energy. Moreover we know that the Drude weight of crs{H = 0) 
has been found to be finite from an approach independent of ours. 



In generic systems the Drude weight, Z), of a conductivity at finite temperatures will 
be zero. It is then the integrability of the 0(3) NLSM and the attendant existence of an 
infinite number of conserved quantities that leads to a finite weight. The existence of these 



quantities can be directly related to a finite D. As discussed in D is bounded from 
below via an inequality developed by Mazur: 

where J is the relevant current operator, Qn are a set of orthogonal conserved quantities, i.e. 
(QnQm) = Snm{Qn)y ^ somc coustaut. For a finite Drude weight we then require that 
at least one matrix element, (JQn), does not vanish. While we do no direct computations, we 
can obtain an indication of whether the matrix elements vanish by examining the symmetries 
of the model. Under the discrete {Z2) symmetries of the 0(3) NLSM, the spin current, J, 
transforms via 

Z^iJ) ^ ±J. 

In order that the matrix element, (JQn), not vanish we require that 

The Z2 symmetries in the 0(3) NLSM include — > —Ua, a = 1,2,3, parity, and time 
reversal. The spin current we are interested in transforms under rotations as a vector. Thus 
any charge, Qn, coupling to the current must also transform as such. From the work by 
Liischer |2^, it is clear there is at least one conserved vectorial quantity such that (JQn) 



does not vanish due to the action of one of the above Z2 symmetries. 
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While the structure of the conserved quantities in the 0(3) NLSM seem to be consistent 
with the existence of a finite Drude weight, this is not the case in the semi-classical approach. 
The dynamics of the semi-classical approximation used in |^,|^ also admit an infinite number 
of conserved quantities (but importantly, different than those appearing in the fully quantum 
model). However as shown there, the structure of the Z2 symmetries in the semi-classical 
approach is such that all matrix elements, (J(5„), vanish. It would thus seem the absence of 
a Drude weight in the semi-classical case is a consequence of differences in the symmetries 
between the semi-classical and fully quantum models. 

To understand the discrepancies in the case of the NMR relaxation rate, 1/Ti, is not as 
simple. However if we believe that the spin conductance demonstrates finite temperature 
ballistic behaviour, it is hardly surprising to find the NMR relaxation rate characterized 
by ballistic logarithms. Again the difference between the fully quantum treatment and the 
semi-classical approach will lie in the differences between the models' conserved quantities. 
Nonetheless one possibility that we must consider is that merely going to 0(e~^^^) in the 
computation of 1/Ti is insufficient. It is possible that we need to perform some resummation 
of contributions from all orders to see the desired singular behaviour, 1/Ti ~ 1/\/H. While 
this would belie our experience with computing the susceptibility and the spin conductance 
via the correlators, the data that goes into the two computations is not exactly identical. 
Thus the possibility that the low temperature expansion of 1/Ti is not well controlled cannot 
be entirely ruled out. 

The differences in the nature of the conserved quantities between the 0(3) NLSM and 
the semi-classical model of suggest the latter is not equivalent to the 0(3) NLSM, even 
at low energies. An indication of this lack of equivalency may lie in the universal nature 
of the ultra low energy S-matrix. This quantity is the primary input of the semi-classical 
model. The semi-classical model imagines a set of classical spins interacting via 

^"^,{9 = 0) = -5ad5cfe, 

i.e. in the scattering of two spins, the spins exchange their quantum numbers. However this 
specification may be insufficient to adequately describe the 0(3) NLSM. Even beyond the 
quantum interference effects which are neglected by the semi-classical treatment, it is not 
clear that the zero-momentum S-matrix is enough to determine the model. 

In this light it is instructive to consider the sine-Gordon model in its repulsive regime. 
The sine-Gordon model is given by the action, 

S = ^j dxdt (^d^.^d^'^ + A cos(/3$) j , (2.24) 

where j3 = P/^/Att. The model is generically gapped. Its repulsive regime occurs in the 
range, in < P"^ < Svr. The model's spectrum then consists solely of a doublet of solitons 
carrying U(l) charge. It is repulsive in the sense that the solitons have no bound states. 
The sine-Gordon model has a similar low energy S-matrix to the 0(3) NLSM, 

Sabi^^ = 0) = -Sad^cb, 

where here the particle indices range over ±, the two solitons in the theory. Thus we might 
expect that sine-Gordon model to possess identical low energy behaviour over its entire 
repulsive regime. 
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This is likely to be in general untrue. For example we might consider the behaviour of 
the single particle spectral function. We might thus want to compute a correlator of the 
form 

(V^+(x,t)v^_(o,o)), 

where ip^ are Mandelstam fermions given by 

/il- i 1 - \ 

^±(x,t) = exp (^± -(-^ +/5)0L(x,t) T -{j- P)^R{x,t)j; 

HiR = \ [^{x, t) ± I £^ dydMy, t)) . (2.25) 

As these fields depend explicitly upon /3, it is hard to see how the properties of the above 
correlator, even at low energies could be independent of this same quantity. More generally, 
(3 determines the compactification radius of the boson in the model and so is related in a 
fundamental way to the model's properties. 

It is useful to point out that Mandelstam fermions are the unique fields that cre- 
ate/destroy solitons that carry Lorentz spin 1/2, i.e. a spin that is independent of 13. They 
would then be the only fields with a chance of matching any semi-classical computation. 
However there are other soliton creation fields, for example, 

for which one could determine the corresponding spectral density. As these fields carry spin 
that varies as a function of j3, their spectral functions will depend upon more than the ultra 
low energy soliton S-matrix. In general, the semi-classical treatment of the sine-Gordon 
model cannot capture its full quantum field content. 

As with the 0(3) NLSM, the conductance of the fully quantum model differs from that of 
the semi-classical treatment. If one were to compute the conductance at finite temperature 
in the sine-Gordon model one would again find a finite Drude weight, D, while the semi- 
classical approach yields D = [^. The notion of under-specificity appears here again. 



The semi-classical approach for the sine-Gordon model equally well describes the Hubbard 
model at half- filling (the solitons are replaced by particle/hole excitations in the half- filled 
band). But it fails to give the correct Drude weight. An analysis of finite size corrections 
to the free energy in the presence of an Aharonov-Bohm flux ||3l| again finds a finite Drude 
weight in the half-filled Hubbard model at finite temperature. 

Interestingly however, there are certain properties at low energies that seem to be in- 
dependent of 13. For example, if one were to compute the low temperature static charge 
susceptibility, the term of 0{e~^^) would be independent of (3. However at the next order, 
C(e~^^^), this would almost certainly cease to be true. And the energy /temperature ranges 
we are interested in exploring do not permit dropping terms of (9(e~^^^). 

It is important to stress we do not question the agreement between the semi-classical 
model and experiment. What we do question is whether the fully quantum 0(3) NLSM ex- 
hibits spin diffusivity. If we are then to understand spin diffusion in terms of the 0(3) NLSM, 
it is possible we need to include additional physics such as an easy axis spin anisotropy 
(weakly present in the experimental system, AgV P2SQ) , inter-chain couplings, or a spin- 
phonon coupling (as done in |TT|). 
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Beyond these, another mechanism that might lead to diffusive behaviour are small in- 
tegrable breaking perturbations of the 0(3) NLSM. Generically any physical realization of 
a spin chain will possess such perturbations, even if arbitrarily small. Such perturbations 
may introduce the necessary ergodicity into the system, ergodicity that is absent in the 
integrable model because of the presence of non-trivial conserved charges, and so lead to 



diffusive behaviour. As discussed in the semi-classical context by Garst and Rosch ^0 
such perturbations introduce an additional time scale, T, governing the decay of conserved 
quantities in the problem. For times, t < T, the behaviour of the system is ballistic and the 
original conserved quantities do not decay. For times, t > T, the behaviour is then diffusive. 
Consequently the Drude weight in the purely integrable model is transformed into a peak 
in cr(co') at C(j ~ 1/T. 

Now the difference in the physics between the 0(3) NLSM and its semi-classical variant is 
not that of integrable breaking perturbations. As demonstrated in their semi-classical 
model is classically integrable. However as discussed above the models do possess different 
conserved charges. It might then seem for certain transport quantities, the semi-classical 
model cures the lack of ergodicity present in its quantum counterpart. 
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III. COMPUTATION OF FINITE TEMPERATURE CORRELATORS 



Here we present the general method by which we compute the correlators at low but 
finite temperature and field: form factor expansions. In the first part of this section we 
consider the general form of these expansions and why we expect them to be applicable at 
finite temperature. In the latter parts of this section, we review the exact expressions for 
the form factors in the 0(3) sigma model together with the necessary regulation of said form 
factors at finite temperature. 



A. General Methodology 

To compute two-point correlation functions, we employ a form factor expansion. At 
finite temperature, such correlators take the form 

G'^{x,t) = ^Tr(e-^^O(x,t)C(0,0)) 

^ En.„ e-/^^-»(n, s^\0{x, 00(0, Q)|n, s^) 

Here t can be real or imaginary time and the sum Ylmsn ^ possible eigenstates of the 

Hamiltonian. Each eigenstate is characterized by the number of particles, n, in the state 
together with a set of internal quantum numbers, in this case the value of Sz carried by 
each particle. The form factor representation of the correlator is then arrived at by inserting 
a resolution of the identity between the two fields: 

e-^^»"(n,s„|O(x,0|m,Sm)(m,Sm|O(0,0)|n,s„) 

At zero temperature, the representation of reduces to one involving a single sum, Em,s,„- 
Thus the computation of amounts to the evaluation of a set of matrix elements. These 
matrix elements can be computed in principle for arbitrary n, m from a knowledge of the 
two-body S-matrix together with various constraints coming from the analytic dependence 
of the matrix elements upon energy- momentum. However with increasing n and m the 
evaluation of these matrix elements and the corresponding evaluation of the sums, En,s„) 
becomes increasingly arduous. 

We are however in a better position when we consider the spectral function corresponding 
to G^: 

E e-^^»»27r5(a; - E,^ + £;,„)(n, Sr\0{x, 0)|m, s^)(m, s^|C»(0, 0)|n, s„) 

= , ^ . 

En.„ e-P^-n (n, s„|n, s„) 

(3.3) 

Wc sec then that only certain terms, those meeting the matching condition, to — Eg^ — Eg^ , 
contribute to the spectral function. 
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In this paper we are concerned in particular with massive or gapped theories. Gapped 
theories are particularly amenable to this sort of computation as they admit a notion of 
thresholds. First imagine fixing n, s„ in the sum above. In a massive theory the intermediate 
states have a finite energy. In particular in the 0(3) sigma model, the energy of an m-particle 
state has a minimum threshold of rriA. And so states with Eg^ exceeding uj + Eg^ do not 
contribute to the sum. For example if uj + Es„ is below the three particle threshold, 3 A, 
states with m > 3 do not make a contribution. 

At zero temperature, i.e. E^^ = 0, the notion of thresholds leads to a situation where 
only a finite number of matrix elements needs to be computed in order to obtain an exact 
result at a given energy, uj. In contrast, at finite temperature we in general would need to 
compute an infinite number of matrix elements in order to arrive at an exact result. However 
here the massiveness of the theory again comes to our aid. With increasing n, the terms are 
weighted with the Boltzmann factor, e~'^^''" < e"^""^. Thus at temperatures small relative 
to the gap. A, we expect in general only the first terms to make a significant contribution. 
We can thus evaluate the correlator in a controlled fashion, expanding it as the sum, 

G^(x,^)=5:c„(x,cu)e-"^^. 

n 

Moreover while the evaluation of this sum in its entirety would require the evaluation of 
an infinite number of matrix elements, each individual coefficient, c„, depends only upon a 
finite number of matrix elements (at least in the cases considered in this paper). As such 
we are able to compute these coefficients exactly. 

While the ability to do so results from each c„ being determined by a small, finite num- 
ber of matrix elements, this feature will not be found in all theories. However form factor 



expansions in massive theories have in general found to be strongly convergent [p!^,[L5Hl7 
Specifically, matrix elements, (n, s„|(9(0, 0)|m, Sm), where n and m are large have been found 
to be relatively small. Even in massless theories where there are no explicit thresholds, con- 
vergence is good provided the engineering dimension of the operator O matches its anoma- 
lous dimension. Thus even if each c„ were determined by a large (even infinite) number of 
matrix elements it would be possible nonetheless to arrive at a reasonable approximation 
for the coefficient. 

There are, however, certain situations where we do not expect to be able to truncate the 
sum, J2n,s„ e~^'^''"( ). In certain circumstances, a physical quantity will see a transition as 
the limit of zero temperature is taken that is non-analytic in nature. To be concrete consider 
the single particle spectral function of the staggered component of the spin field: 

S{x,u) = {n{x,u)n{0,0)). (3.4) 

At zero temperature, we expect that for energies, u < 3A, S{x,uj) takes the form of a 
5-function: 

S{x, u) = cS{u - A), u <3A, T = 0. (3.5) 

However at finite temperatures this (5-function is broadened. We then do not expect to be 
able to see this broadening unless we evaluate the sum, J2n,s^ e"^'^"" , in its entirety. Indeed, 
computing S{x, u) through the truncation of this sum at any finite n leads to 
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S{x,u) = cS{uj-A) + --- T>0. 



(3.6) 



Only through the resummation of the higher order terms is the 5-function replaced by a 
broadened peak. However it may well be possible to guess at the resummation on the basis 
of the first terms in the series. 

Rather than consider such situations, we want to focus upon quantities that possess a 
smooth T transition. As such consider the behaviour of the staggered field spectral 
function, S{x,uj), for energies below the gap u < A. At T = we have 

S{x, \uj\ < A) = 0, (3.7) 

while at T 7^ 

S{x, \u\ < A) = C(e-^^). (3.8) 

Thus the T — > limit behaves in a smooth fashion. 

This method is markedly different than that developed in P3|-p7[. In our method we 



employ the basis of eigenstates, \n, that arises from the zero temperature problem. There 
a new basis is adopted that takes into direct account the thermalization of the vacuum state. 
Let I Or) be the state with a representation of the particle content of the system in equilibrium 
at finite T and let |(n, s„)r) be states that are excitations above this thermalized ground 
state. (In contrast, \n, s„) are excitations above the empty vacuum state.) With such a 
basis, the correlators have the following form factor representation: 

(O(x,t)O(0,0)) = 5](0T|O(x,t)|(n,s„)T)((n,s„)T|O(0,0)|0T). (3.9) 

n,s„ 

This method involves considerable technical complications. In general, it is a challenge 
to compute the new vacuum state \0t) as well as the excitations above \0t), never mind 
the form factors (Ot|(9(x, t)|(n, s„)t)- These difficulties are only enhanced by the non- 
diagonal scattering present in the 0(3) sigma model, i.e. the two-body S-matrix is other 
than S^jj'^ = Saa'^bb', where no internal quantum number are exchanged. This method was 
developed in particular for theories that are massless. However in our case the theory is 
gapped. It thus makes sense to exploit the control over the sum, J2n,sn e'^^"^ , that the low 
temperature regime affords us. 



In some sense our approach is similar to that of LeClair and Mussardo [IS]. There they 
begin with the form factor sum as in ( |3.1D . However they recast the sum of the thermal 
trace through introducing a set of hole excitations complementary to the particles. Hole 
excitations appear naturally in terms of the form factors. A typical form factor that needs 
to be evaluated for a finite T correlator looks as follows, 

(si,e,JO(a;,t)|s2,e,,), (3.10) 

where we have explicitly labelled the energy of the state. Using crossing symmetry, this 
matrix element can be rewritten as 

(0(a;,t)|s2,es2; si, -e^,), (3.11) 
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provided ei = e2, Si = S2 does not hold. Here Si is the 'charge conjugate' of si. The 
excitation, (si,— i^^J, can be thought of as a new type of excitation, a hole. Thus the 
double sum of a two point correlator was recast in |]19| as 



(C»(x, t)O(0, 0)) = J2 n fK) n fi^sj{0{x, t) Imp, Sp- m^, s;,) (m^,, s;,; m^, Sp|O(0, 0)). 



(3.12) 



Notice that the partition function, Z, is absent from ( ^.121 ) while new factors, Ylf, have 
been added to the expression. Each /(e^) is the occupation number of the excitation (in this 
case assumed to be fermionic), s, with energy e: 



These modifications represent an ansatz put forward in 
the regulation of the matrix elements, 

(si,ei|0(x,t)|s2,e2). 



T^, and are argued to come from 



(3.13) 



in the case si = S2, ei = 62- 

Although this ansatz is supported in the case of one point functions (i.e. expectation 
values of the energy or spin) p9|-pT| , it has come under criticism for the computation of two- 
point functions in [^. There the allied case of current-current correlators in the quantum 
Hall edge problem at T = but finite voltage was examined and it was found that their 
ansatz did not seem to reproduce the correct results. 

What relevance does this critique have for our approach? We do not and cannot use the 
ansatz of LeClair and Mussardo as scattering in our theory is non-diagonal and their ansatz 
only makes sense in the case of theories that are diagonal. However might the critique in 
pO| still have bearing upon our results? We do not think so. The correlator considered in 



20| is computed in a massless theory whereas our results depend upon the gapped nature 



of the 0(3) sigma model producing a series of thresholds. Moreover we already expect to 
run into difficulties whenever there is non-analytic behaviour at T = near a threshold as 
in the behaviour of the staggered field spectral function near a; ~ A. Thus we do not expect 
to capture the physics of the conflation of all the thresholds in a massless theory. 



B. Form Factors in the 0(3) Sigma Model 

1. Constraints Upon Form Factors 

The form factors of a field O are defined as the matrix elements of the field with some 
number of particles, Aa{0): 

fZ-aSOu ■■■,0n) = {O{0, 0)A^M ■■■Aa, (^l))- (3.14) 

The Aa{6) are Faddeev-Zamolodchikov operators which create and destroy the elementary 
excitations of the theory. 6 is the rapidity which encodes the energy-momentum carried by 
the excitation. 
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p=Asinh(6'); £" = A cosh(6'). 



(3.15) 



The form of fai---an determined by a combination of two-body scattering, Lorentz invari- 
ance, analyticity, and hermiticity. 

The constraint from scattering is derived from the commutation relations of Faddeev- 
Zamolodchikov operators: 

= szti(^^-e2)Al{eMliOsy, 
AlieMaM) = KaAdi - O2) + szz^Oi - e,)A,,{e,)Al{e,). (3.16) 

S, the two-body S-matrix, gives the amphtude of the process by which particles {01,02} 
scatter into {03,04}. It is solely a function of 61 — 62 = 612 by Lorentz invariance. In our 
case, scattering between magnons in the 0(3) model, the S-matrix is given by 

2ni9 



cr2i9) 



+ in){e - z27r)' 
9{e-iTT) 



_ 27ii{in - 9) 
[0 + m)[6 — zzTT) 

As 9 ^ 0, the S-matrix reduces to S^^^t ~ ~^aia4,Sa2a3- This is the approximation underlying 
the semi-classical analysis of Damle and Sachdev J^,^] . For the form factor to be consistent 
with two body scattering we must then have 

Si;;:^'iO. - 0.+,)fZ...,a'.,a'.,,,...,aAOu ■■■,9,, 9,+„ ■■■,9^). (3.18) 

This relation is arrived at by commuting the i-th and i + 1-th particle. 

A second constraint upon the form factor can be thought of as a periodicity axiom. 
In continuing the rapidity, 9, of a particle to — 27ri, the particle's energy-momentum is 
unchanged. However the form-factor is not so invariant. We instead have 

fZ-,a. (SU ■■■,0n)= fZ,a,,-,a„.. i^n - 2m, ^i, ■ ■ ■ , 9^.,). (3.19) 

This constraint is derived from crossing symmetry [113. It implicitly assumes that the field O 



is local: if O is non-local additional braiding phases appear in the above relation p0| , |39| , |35 

Another condition related to analyticity that a form factor must satisfy is the annihilation 
pole axiom. This condition arises in form factors involving a particle and its anti-particle. 
Under the appropriate analytical continuation, such a combination of particles are able to 
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annihilate one another. As such this condition relates form factors with n particles to those 
with n — 2 particles. In the case of the 0(3) sigma model it takes the form 



i reSe„=e„_i+7ri/(6'l, ■ ■ ■ , dn)a^,-,ar, = f{di, ■ ■ ■ , 9n-2)a[,-,a'„_2^a^a'^_-^ 

— Sriai ^ {On- ll)Sll'^^ {9 n- 12) " " " 5'T„_3a"_|(6'.„_i„_3)S'a„_ia^_2 (^ri-ln-2)y (3.20) 

This relation as written assumes that we are normalizing our particle states as {0\6') = 
2n5{e-e'). 

The form factor must also satisfy constraints coming from Lorentz invariance. In general, 
the form factor of a field, O, carrying Lorentz spin, s, must transform under a Lorentz boost, 
9i ^ 9i + a, via 

C-a„ (^1 + «,■■■, e„ + «) = e^VZ-a,. (Oi,---, e„). (3.21) 

The particular fields we will be interested in are the magnetization density, Mo{x, t), as well 
its corresponding conserved current, Mi{x,t). Together they form a Lorentz two-current. 
(Here 0, 1 are Lorentz indices. Spin indices have been suppressed.) As this current is 
topological we may rewrite it in terms of a Lorentz scalar field, m{x,t): 

M^{x,t) = e^,d''m{x,t) (3.22) 



The form factors are then determined for the field m{x,t) which obeys ( 3.21|) with s = 



while the corresponding form factors of Mf^{x,t) are related to those of m{x,t) by 

/a^^J^i, ■■■,0n)= e,.P-'{e,)fZ...aAOu ■ ■ ■ , ^.), (3-23) 



where P° = Acosh(^i) and = Ei Asinh(^i). 

These conditions do not uniquely specify the form factors. It is easily seen that if 
f{9i, ■ ■ ■ , On)ai,-,a„ satisfies these axioms then so does 

f(n n N Pn(cOsh(%)) 

Q„(cosh(%)) 

where P„ and Qn are symmetric polynomials in cosh(^jj), 1 < i, j < n, and are such that 

Pn\e„=e„^i+TTi = Pn-2] Qn | e„=0,j_i+7ri = Qn-2- (3.25) 

To deal with this ambiguity, we employ a minimalist axiom. We choose P„ and Qn such 
that Pn/Qn has the minimal number of poles and zeros in the physical strip, Re(^) = 0, 
< Im^^ < 2-71. Additional poles are only added in accordance with the theory's bound state 
structure, an unnecessary complication in our case as the 0(3) sigma model has no bound 
states. Using this minimalist ansatz, one can determine Pn/Qn up to a constant. 
To determine this constant we rely upon the action of the conserved charge 

S^= [ dxM^{x,0), 
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upon the single particle states. We expect 



{9, s, = i\s,\e\ s, = i) = 2Tr6{e - e'). (3.26) 

Thus from the knowledge of the two particle form factor, we can fix the overall normalization. 
To check this normalization we will compare the form factor computations with the results 
of other techniques. For example we will compute the magnetic susceptibility using both 
form factors and the thermodynamic Bethe ansatz. Through comparing the T — > 0,H —>■ 
results, we see that the normalization has indeed been consistently computed. 

As yet another check we can fix the phase of this constant using hermiticity. For this 
purpose it is sufficient to consider 2-particle form factors. Hermiticity then gives us 

{0{o,o)AaM)AaMr = «(^iX(^2)ono,o)) 

= {O\0,0)Aa-,{9i - in)Aa-M - zvr)), (3.27) 
where the last line follows from crossing and so 

C(^i,^2)* = - in,e^-tn). (3.28) 



C. Review of 0(3) Sigma Model Form Factors 



From ( ^.22 ) and (^.23) it is sufficient to give the form factors for the scalar operator. 



m{x,t). These have been computed by both Smirnov [|T4| and Balog and Niedermaier p2 



However |3^ presents them in a more amenable form, possible in this particular case because 
of the simple structure of the S-matrix of the 0(3) sigma model. 



Using the axioms as presented in the previous section, thus finds for the two and 
four particle form factors 



121314 ' 



nV'(%) X (5''*''^e""^"^^i(^i) + (5''*''^e""^''ic/2( 

i<j 



^^a4m^aasa2g^(^Q,^ + 6''''''^ e"''^''' gi{9,) + d"'-""' e""^"^ Q^iOi) + 5'^2"ie"''^"«^6(^*) 



V^6(^, 



(6^32 - i7r)6'3i(6'3i - m) 
(^32-^vr)(^32 + ^27r)(m-03i) 

^32^3l(3vr2-03l) 
6*32 (6'32 - ^7r)6'3i 

2'Ki{m - 632)031 
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+i(6'43 - «vr) 



m{6^2 + 6*31 
-27r2 - Syria 



31 



+ 9: 



732 
2 

31 



^31)' \ 



-47r2+m(a32-2^^3l)-^32' 
2tX^ + Z7r(^^32 + 2^31) - 2^^32^^31 
-^7r(2a32 + ^31) + 2^^32^^31 



'43, 









732 



^31 - 27ri 

\ ^32 — ^31 / 



(3.29) 



We have checked that these form factors do indeed satisfy the necessary axioms and found 
that the results of |^ are without typographical error. The reader should note however 
that we use a different particle normalization than |^2|] and so the results differ by an overall 
multiplicative constant. 

The two particle form factor differs from that appearing in Affleck and Weston's work 
^ on the 0(3) sigma model. The two particle form factor Affleck and Weston use is given 
by 

xA/f rn n\ ( ^ ln\ ( n \\ 3ai 09 ^^^^(^^12/2) ZTT + ^^12 
faX(^l,92) OC (cosh(ai) - C0sh(a2))e'"'"' ' ' -^n " 

This differs from our form in that it has a different Lorentz structure and lacks an extra 
factor of tanh(6'i2/2). The different structure of the two particle form factor is a result of the 
constraint the annihilation pole axiom places on form factors of different particle numbers. 



If one only computes the two particle form factor, as done in [Q, this constraint can go 
unsatisfied. However in terms of the low energy behaviour (i.e. 61,62 —>■ 0), the two forms 
for f^a2 nearly identical. 



D. Regularization of Form Factors 

We end this section with a discussion of the regularization of form factors that appear 
in the evaluation of thermal correlators. Form factors with all particles either to the right 
or the left of the field such as 

do not pose any such problems. However the form factors encountered in the evaluation of 
finite temperature correlators are of the form 

To understand such an object we must contend with the possibility that 6i = 6j, ai = bj 
for some From the algebra of the Fadeev-Zamolodchikov operators ( |3.16| ), we know the 
commutation relations involve 5-functions, i.e. 



Al{6,)A,{6,) = 27c6{6,-6,)Sa^,^ + 



(3.30) 
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It is crucial to include the contributions of the (5-functions to the correlators. In particular 
they contribute pieces which cancel off otherwise ill-defined terms arising from the partition 
function. To do so we must understand the above form factor to equal 

{AbjL) ■ ■■A,,ie,)Oio,o)Aa^ie^) ■ • ■ ^,,(^0) = 

E SA,AiSB,BABl\Ai){B2\OiO,0)\A2)connected. (3.31) 

{aJ=AiUA2 

The sum in the above is over all possible subsets of {ui} and {bi}. The S-matrix Sa,Ai arises 
from the commutations necessary to rewrite Aa^{9n) ■ ■ -Aq^ (6*1)10) as A2A1IO) and similarly 
for Sb,Bi- The matrix element {Bi\Ai) is evaluated using the Fadeev-Zamolodchikov algebra. 
In this way (ill-defined) terms proportional to 5(0) are produced but which cancel similarly 
ill-defined terms arising from the evaluation of the partition function. 

The 'connected' form factor appearing in the above expression is to be understood as 
follows. Using crossing symmetry, the form factor can be rewritten as 

(52|O(0,0)|A2)connectcd = (A' (k) ' ' ' ^b' {e,,)0 {0, 0) A,>, (% J " " " ^a' (% ))connected 
= (0(0, 0)Aa' (^jj ■ ■ ■ Aa' {ej,)Ai, {9,^ -in)--- Ay_ {Q^^ - m)) connected 
~ fb'. ■■■b' a' ■■■a' i^h ~ '^'^ ■• ' ' ' ' ^«fe "^^5 ^ii ' ' ' ' ' ) connected , (3.32) 

where the last relation holds provided we do not have 9i = 6j, ai = bj for any If this 
does occur we see from the annihilation pole axiom that the form factor is not well defined, 
having a pole at 9i = 9j. In such cases the form factor requires regulation. 

To regulate the form factor, we employ a scheme suggested by Balog P3 and used by 



LeClair and Mussardo 191. We define 



f'b' -V a' -a' i^h - ^71" + ^^71, ■ ■ ■ ,9i^^ - in + ir]k, 9j^,--- , 6*^) connected 
'1 'k h Jq 

= finite piece of lim ...5, , . , {9i^ - in + ir]i, ■ ■ ■ ,9i^ - in + ir]k, 9j^,---, 9j). (3.33) 

In taking the finite piece of /*^, we discard terms proportional to r]~^ as well as terms 
proportional to rji/rjj. In this way the connected piece is independent of the way the various 



limits rji are taken. Balog |^3[ has already used this prescription to compute one point 
functions and successfully compare them to TBA calculations. In [^] it was argued that the 
delta functions leading to such terms arise from the use of infinite volume wavef unctions. If 
such wavefunctions are replaced instead with finite volume counterparts, the delta functions 
are regulated. For example, a pole in r] is changed as follows 

l^l^e^^l^m ,3.34) 

where f{9) is some sharply peaked function about ^ = which in the infinite volume 
limit evolves into a 5-function. However the principal value of this regularized integral is 
zero. Thus discarding the pole terms is justified in this sense. For terms that are ratios of 
infinitesimals, Balog also demonstrates that such terms, once regularized, disappear in the 
infinite volume limit. 
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IV. THERMODYNAMIC BETHE ANSATZ AT FINITE TEMPERATURE AND 

FINITE FIELD 



In this section we review the derivation of the equations describing the exact free energy 
(and hence the susceptibihty) of the 0(3) sigma model together with its low temperature 
expansion. The exact description of the thermodynamics of the 0(3) sigma model takes 
the form of a set of quantization conditions for the momenta, Pa, of the excitations in the 
ground state. With = Asinh(^cK), we have the following condition 



A. 



N n 

TT- 



da — Xa + «vr 



N 

^1 Op - Xa- IT^ 



-Q A^ — Aa + ivr 



7=1 -^7 



A. 



(4.1) 



Here Op are the rapidities of the other excitations in the ground state while the A's mark out 
spin excitations above an originally polarized ground state. (The Bethe ansatz construction 
begins with a completely polarized ground state of spin 1 excitations above which one then 
creates spin excitations - marked out by the A's - in order to give the ground state the 
desired spin polarization.) is the total number of excitations in the ground state while M 
is the number of spin excitations. The quantum number, Sz-, is then given hj Sz = N — M. 
The total energy of the ground state in a magnetic field is then equal to 

N 

E = AJ2 cosh(^a) - H{N - M). 

a=l 

The analysis of these equations proceeds using the string hypothesis. The solutions of the 
above equations take the form 

On are real; 



A 



n,k 



Xl + i7i{n + l-2k)/2, k 



1,2, 



(4.2) 



that is the A's are organized into 'complexes' which share a real part, A^, the centre of the 
complex. 

In computing the free energy, we are interested in the continuum limit of the above 
equations. To arrive at this limit, we introduce densities per unit length, p{0) and cr„(A), of 
respectively the O^s, and the centers. A", of the complexes. We further introduce particle 
and hole densities by writing p = + pp and (j„ = anh + o'np- A particle density gives the 
probability that the ground state contains an excitation at a given rapidity, 0/X, while the 
hole density gives the converse probability that the excitation at the rapidity is not found 
in the ground state. Equations describing these densities can be arrived at in a standard 
fashion (see Section 8.3 of for an analogous derivation in the case of the Anderson 
model): 

pp{0) + ph{0) = ^ cosh(e) + [s * a^hM; 
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0"mp(A) + (TmhW = S2mS * Pp(A) + S* {(Tm+l,h + 0-m-l,/i) (A) , (4.3) 

where s{x) = (vr cosh(x))~^ and f * g denotes the convolution of these two functions: 

f*g= f dX'f{X-X')g{X'). 



From these equations the free energy per unit length, can be derived (again see |36] for 
details of an analogous derivation): 



TA 



d9 cosh(9) log(l + e 



e{e) = Acosh(^) -Ts*log(l + e^'2)(^); 



e„(A) = Ts * log(l + e^^"-i)(l + e^^"+^)(A) + 62nTs * log(l + e~^^)(A); 



lim — = H. 

n — >-oo fl 



(4.4) 



Here we have expressed the free energy of the system in terms of the dressed energies (or 
pseudo-energies), e/e„, of the excitations. These functions give the energetic cost of making 
an excitation at a given rapidity taking into account the excitation's interactions with the 
other particles in the ground state. These equations are in agreement with where they 



were first written down and correct typos found in [11 



To derive the low temperature expansion of the free energy, we follow ||33| • We solve the 



above equations 0| through iteration. We write for each pseudo-energy, 



1 + e 



(4.5) 



m=0 



This expansion is such that r„m is of 0{e "^/^^). On the basis of (|4.5|) we can write the free 
energy as a series in e"™^^: 



(4.6) 



m=l 



We will compute the m = 1, 2 terms of this expansion. 

The m = term of ([4.5| ) is arrived at by neglecting the term involving log(l + e"'^'^) in 
the equation for e„. If this is done, these equations reduce to 



T 



e„(A) = ;-log(l + e^^"-i)(l + e 



lim — = H. 

n— >oo fl 



(4.7) 



They are then algebraic in nature and admit the following solution: 
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sinh(f (n + 1)) 

At this order of the iteration, e{6) becomes 

e{e) = A cosh(0) - T log 0(2) , (4.9) 

and so 

TA r 

n = -— / d^cosh(^)log(l + 3e-'^^=°^^(^)). (4.10) 

ZTT J 

Clearly Q is of C(e-^^). 

The next coefficient in the series ( |4.5| ), r„i, is found by substituting (O) into the equa- 
tions for e„: 

e„ = Ts * log(l + e^^"-^)(l + e^'-^') + 62ns * log(l + 0(2)e-'5^™^^^(^)); 
lim — = if. (4.11) 

n— >oo 

To the order in e~^'^ to which we are working, these equations reduce to 

"^'^"^ Tni = (rn-1,1 + Wi,i) * s + 52nTlog(l + e-^^^°^^(^)) * s. (4.12) 



(/)(n - l)0(n + 1) 
As can be directly checked, they admit the solution 
0(1) 



r-ni 



(0(n + l)a„ - 0(n - l)a„+2) * * T log(l + 0(2)e-^^^°^^(^)); 



2n 

With this, e{9) to 0(6"^'^) takes the form 
e{e) = Acosh(^^) - riog0(2) - T^(0(3)a2 - ^(1)04) * e-^^^°^^(^) + 0{e-'^^). (4.14) 



We can continue this procedure, obtaining rnm,'^^' > 2. Indeed goes on to compute r„2 
and so corrections of 0(e~^^^) to 6(6*). 
The zero field susceptibility is given by 



;^(// = 0) = -d^H^\H=0 = -7r f C?^COsh(^)-^ 

ZTT J i 



j/=o 



(4.15) 



Using e(^) in ( |4.14|) and expanding the above expression to 0{e ^'^'^), we find 



^(^H = 0) = (^ I (i^cosh(^)e-^'^™^'^(^)(l - 3e-^^ 



cosh(9)^ 



/ rfgirfg2Cosh(gi)e-^-"(-^(^-)+-^-(^^)) ^, ti^'"'. , +0(e-^^-"), (4.16) 



'12 "T ^12 

where 6'i2 = 61 — 62. This agrees exactly with the derivation of x coming from the compu- 
tation of the two and four particle form-factors. 
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APPENDIX A: COMPUTATION OF MAGNETIC SUSCEPTIBILITY USING 

FORM FACTORS 



To compute the correlator, {Mq{x, t)Mq{0, 0)), we first consider the action of the thermal 
trace: 





|M3(x,r)M3(0,0)|n,s„) 




\n, Sn) 



C{x, t) = (Mo^(a;, r)Mo^(0, 0)) = ^ ^ J.^. „ v " ■ (Al) 

Keeping the first two terms leads us to 

(M3(a;,r)M3(0,0)) = (| ^e'^^^-^^W ^ (yl„(^)|Mo^(a;, r)Mo3(0, 0)|A„(^)) 
1 r dOi c?^2 /3A(cosh(ei)+cosh(e2)) ■ 



2 J 27r 27r 



Y,{AaM)Aa,{e2)\M^i^,r)M^iO,0)\AaMA,^{ei)) 
+ ^e-^^-M^)(A„(e)|A„(e))) (A2) 



Expanding the denominator then gives us 

d6 



1 f dOid92^_pA{cosh{0i)+cosh{e2))^ 

2 J 27r 27r 



0102 



The term arising from the partition function is ill-defined as the state normalization is given 

by {Aa{6)\Aaj^{6i)) = 27TSaaiS{9 — Oi). However this term will be canceled by disconnected 

terms arising from {Aa,{e^)Aa,{e2)\Ml{x,T)Ml{QMAaM)AaM)- 

To evaluate this expression we begin by computing the first term of the trace 
(Aa(6')|MQ (x, r)MQ (0, 0)|Aa(6')) by inserting a resolution of the identity between the M^'s: 

dOi f dOn 



x{AM\Ml{x,T)\AaSen) ' • • ^.K)) Ri (^i) " |M^(0, 0) . (A4) 

We only need to keep the lowest order term, n = 1, in this expansion; all other terms make 
no contribution to the susceptibility. The terms corresponding to n even are identically 
zero (by parity); the remaining n > 1 odd terms vanish in the low energy-low momentum 
limit of the corresponding spectral function. (We will return to this is a moment.) Given 
that we can thus compute the entire contribution (y4a(6')|MQ (x, r)MQ (0, 0)|Aa(6')) makes 
to the susceptibility, we will able to find an exact correspondence between the form factor 
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computation and a low temperature expansion of the exact free energy. With the n — 1 
term we then have 

{AM\Ml{x,T)Ml{QMMO)) = T. j ^{AM\Ml{x,r)\A^M)) 

x{A,^{e^)\MlM\A,{e))- 



^^lg-rA(cosh(ei)-cosh(6»))+ia;A(sinh(6»i)-sinh(e)) 



^ [Ml 

x{M^{0,0)\AaMAa{9-in)){M^{Q,0)\Aa{9)AaAOi-i7r)); 

^Y^j ^g-rA(cosh{eO-cosh(e))+ixA(sinh(ei)-sinhW) yM|(^ _ - .^T, 9). (A5) 

We have used crossing symmetry in the second hne. From Section 3.3, the form factor 
fail {9, 9i) is given by 

faaliO, Oi) = i^e=^-Hsmh(^) + sinh(^i))V'(^ - ^i). (A6) 
Then the lowest order contribution, Ci{x, r), to the spin-spin correlator, C{x, r), is given by 



cosh(e)g-rA(cosh(9i)-cosh(6»))+ixA(sinh(6li)-sinh(e)) 

xf^J{9-i7^,9,)fil^{9^-i7^,9). (A7) 



Fourier transforming in x and r and continuing cun — —iuj + S yields, 

Ci{uj = 0, A; = 0) = ^ / c^^cosh(^)e-^^'^°^^(^) = ^Xi(/3A), (A8) 

where Ki is a modified Bessel function. This has the expected small temperature behaviour. 

Let us consider further why the above computation gives the sole contribution to 
(Aa(^)|Mo(x)Mo (0)|Aa(^)). The next potential contribution to this matrix element takes 
the form 

J d9id92d9s{Aa{9)\M^{x)\AaA9l)Aa,{92)AaM)) 

{A,, (^3) (92) A,, {9i)\Ml{0) \A,{9)). (A9) 

Upon evaluation this expression produces two types of terms. The first is associated with 
the disconnected pieces of the matrix elements appearing in the above. An example of this 
type of term is 

= J d9id92e'^''^'^^^^'^+'^^^^^^\smh{9i) + sinh(e2)) x (term regular in 9i and ^2)- (AlO) 

34 



As Mq is a Lorentz current, the term (sinh(^i) + sinh(^2)) appears in the above. Thus 
when the Fourier transform, / e*'^^, is taken followed by the limit, k ^ 0, this term vanishes 
identically. 

The second type of term we must deal with in evaluating (|A9|) takes the form 

J d^id^2rf^3e'"^(""'^('^)+"°'^('^)+"°'^('^)-"°'^(')VfLai(^ - «vr + te, 63, 62, Oi) 

xfaXami^i -i'rr + iei, 62 - in + zes, 63-171 + zeg, 9), (All) 



and arises from the connected pieces of the matrix elements appearing in ([A9|). To evaluate 
this term we deform^ the contours ^1,2,3 via 

^1,2,3 (^1,2,3 + 

In doing so we deform through a number of poles whose residues we will pick up. Evaluating 
these residues we again obtain something of the form ([A10|) . As such, Fourier transforming 
and taking the k = limit forces the term to vanish and ( |A9| ) ends up making no contribution 
to the susceptibility. In the same way, it is easy then to convince oneself that terms involving 
an even greater number of particles similarly do not contribute to the static susceptibility. 

We now go ahead and compute the second term arising from performing the thermal 
trace, {Aa,{6i)Aa^{e2)\M^{x,T)M^{0,0)\Aa2{02)Aa,{6i)) . We evaluate it as before by in- 
serting a resolution of the identity between the two fields. In this case the only term that 
contributes is the n = 2 term: 



{AaA0l)Aa,i92)\M^ix,T)M^{0,0)\Aa,ie2)AaA0l)) = 

E J -^^{AaM)AaM)\M!{x,T)\A^,{e3)AaM)) 



x{A^M)AaM)\M^iO,0)\AaM)AaM)- (A12) 

Allowing for the presence of disconnected terms, the matrix elements in the above expression 
take the form 

{AaMAa,{92)\M^ix,T)\A,M)AaM)) = 

6a,a,2n6{d^-d,)ffS^{e2-tn,e3) 



^In doing so we assume that time is real, not imaginary. This does not pose a problem as we could 
as well directly evaluate the retarded correlators as opposed to evaluating them via an analytical 
continuation of imaginary time-ordered correlators. 
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+ /a2°ai,a4,a3(6'2 - «71", Ol " ^TT, 6*4, 6*3)^ (A13) 



where fc refers to a connected form-factor. We now substitute ( |A13|) into ( A12 ) and obtain 
the following after some lengthy but straightforward algebra 

4=aia2a3aJ StT 271 2tX 2tX 

X {Aa, {e^)Aa, (^3) |M3(0, 0) I (^2) Aa, (^l)) 
= C2I + C22 + C23 + C*24 + C25 + C'26; 



C2i = 27r5(0) 5: I 



27r 27r 27r 



3 g - /3 A (cosh (01 ) +cosh (6I2 ) ) 



^^-TA(cosh(03)-cosh(e2))-tx'A{sinh(03)-sinh(6l2)) 
X c 



--y _ _o I C^^l '^^2 2/3Acosh(6)i) -TA(cosh(6>2)-cosh(ei))-izA(sinh(6>2)-sinh(ei)) 

~ ' 27r 27r 

aia2 

^ _ j '^^1 '^^2 /jA(cosh(ei)+cosh(6>2)) -rA(cosh(6>i)-cosh(6>2))-ix'A(sinh(ei)-sinh(6>2)) 

~ ' 27r 27r 

X E -^^>^2)/£l(^2 -^vr,^i); 

aia2 



^ ^ 1 /■ ^02 '^^3 /3A(cosh(6>i)+cosh(6>2)) 

A„^„J 27r 27r 27r 



111213 



^„-rA(cosh(6»3)-cosh(6»2))-«xA(smh(6l3)-sinh(6»2)) 
X c 



X |/ijf°ai,ai,a2(^3 " ^TT, 6^1 - OT, 6*1, 6'2)c/&3(^2 - ^VT, 6*3) 



+ E ^a^1(^2l)^:2^:K^13)&,a4,a3(^l " ^3 " ^2, ^l)c/,:>^(^2 - ^VT, ^3) 

+ (^2--^^3)|; 



n ^_ [ ^^1 ^^2 ^^3 /3A(cosh(gi)+cosh(6'2)) 

4,^i;'a3^ 27r 27r 27r 

^^-rA{cosh(03)-cosh{6'2))-23:A(sinh(6'3)-sinh(^2)) 

A 



X<| E 5a1:^(^2l)/S,a,,a3(^3 -^Vr,^! -Z7r,02,ei)c4l(^2 -^7r,03) 
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+ E Slttl (^^13)/S.a4,a3 (^1 - ^3 - ^^T, ^^i, ^2)0/^';,, (^2 - ^VT, ^3) 
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^ ^ 1 f ^^^^ -/3A{cosh(ei)+cosh{02)) 



aia2asa4 

_ ^-TA(cosh{6'3)+cosh(6*4)-cosh(6'i)-cosh(02))-za:-A{sinh(6'3)+sinh(04)-sinh{6'i)-sinh^^^ 

X/a2°ai,a4,a3(^2 - ^VT, Oi - in, 6^, ^3)0- (Al4) 

Although appearing exceedingly complicated, these terms dramatically simplify once we 
Fourier transform. 



The first term, C21, on the r.h.s. of ( A14 ) involves 5(0) and so is ill-defined. However it 



precisely cancels the term arising from the evaluation of the partition function in 

1^1^ E e-^^(™^^(^^)+^°^^(^^»(A„,(^)|Mo^(x,r)Mo3(0,0)|A,,(^))(A,,(^ (A15) 

27T 

as is evident if a resolution of the identity is inserted between the two fields, Mq, in the 
above and then truncated at the one-particle level. 

Having canceled off the 5(0)-terms we now look at terms that make a genuine contribution 
to the spin-spin correlator. We first consider the completely disconnected terms. Fourier 
transforming C22 and C23 in time and space and then analytically continuing, Un — > —iu + S, 
leads to 



3 A r°° 

C22{uJ = 0, = 0) + C23(cj = 0, A; = 0) = -3— / rfecosh(^)e 

TT J-00 



-2/3A cosh{e) 



= -6— iri(2/3A), (A16) 
vr 

where again Ki is a standard Bessel function. 

To compute C24 we need to evaluate the connected four-particle form factor. To do so 
we add small imaginary pieces to the rapidities where potential poles lurk and take only the 
finite piece. For example the first term of C24 upon Fourier transforming reduces to 

C2i{uj = 0, A; = 0) = / rf^iC^^2e-^^(™'^^('^)+'°^^('^^) cosh-^(^2) 

^f2ihd2 - 62) E /i^ai2(^^2 - ^VT, 9^ - ITT, 6,, 62), + three other terms. (A17) 

ai 

Then to evaluate the connected form factor in this expression we write 

/laiai2(6'2 " ^TT, 6*1 - ZTT, 6*1, 62)0 = 

finite part of /ia°ai2(^2 — "^ti", — "^tt, 9i — it], 62 — i5) (A18) 

We evaluate this matrix element using the discussion in Section III.D, throwing away any 
poles in ?7 or 5 together with terms of the form ri/5. Expanding the form factor on the r.h.s. 
of ( [A18| ) in T] and 5 by using ( |3.29| ) leads to 
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fm°aA^2 - «7r, 9i - in, Oi - ir], 62 - i6) 



ai 



Avr^ 16 COsh(0l) COSh(^2) w-pr , , .^m^ in - a ■ n n\ 
15 )(ll^l5=0,„=o)G'laiai2(^2 " - 27r, , ^2) 

+ coMd2)d^r,[(RmTlaA(^2-i'n.e,-in,e,-i^,e,-i5)Y^^^^^^^ (ai9) 

where H 'V' is given in this case by 

Jl ^ = ^(^2i)V'(^2i - ivr + iri)ip{ei2 - ivr + i5)^(^i2 - ir/ + 

Discarding the pole terms (the first set of terms on the r.h.s. of ( [A19| )) and evaluating the 
remainder leaves us with the desired connected form factor 

Combining ([K20D and (p^ with (^T^) we find 

^24(0; = 0, = 0) = ^ / rf0irf^2e-''^(™^^(''^)+^°^"('^» 

6n^ cosh(^i) + (Stt^ + ) coshf^a) 

" i4.^ + eh){n^ + eh) + 

11 A/? T 
= —4-Ko(pA)Ki(pA) + 0(-e-f^^) + three other terms. (A21) 
An-^ A 

To arrive at the last line we have dropped terms polynomial in 612. This leads to errors of 
0{^e~^^). The remaining three terms make equal contributions to C24. We thus finally 
have 

^24(0; = 0, = 0) = ^ / ^0^^0^^-mcosHei)+cosHe,)) 



^^ 67r^cosh(gi) + (57r^ + 2g^2)cosh(g2) ..... 

(An^ + ei2)in^ + eh) • ^""''^ 

We note that in regulating the form factor for C24 we do not allow the infinitesimal imaginary 
pieces to affect the spatial dependence of the form factor, i.e. we do not write 

If we were to do so we would find an additional term coming from expanding exp(iAx • • •) 
in 7] and 5. However generically such terms lead to a violation of translation invariance 
and as such should not be included. We moreover know that such terms would violate the 
equivalence of the form factor computation with the expression for the susceptibility coming 
the thermodynamic Bethe ansatz. 
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We go through a similar procedure with C25 and find an identical result: 025(00 = 
0,k = 0) = C2a{ijJ = 0,k = 0). That we do so is significant. We might have approached 
the calculation equally validly by ordering the in and out states in the thermal trace and 
resolution of identity such that 6*1 < 6*2 and 6*4 < 9^ (and correspondingly multiplying the 
expressions in ( [A14| ) by 4). If we had done so we would find that in this case C25 = and 
C24 is twice its current value. Of course both approaches must yield the same answer. But 
to do so we need 6*25(0; = 0,k = 0) = 6*24(0; = 0, = 0). Given the regularization of the 
form factors one must do to compute 6*25, it is not a priori that this will be the case. That 
it is is a non-trivial check of our regularization procedure. 

The final term we must evaluate is C*26. Fourier transforming as before we find 

1 f dOi (162 de^ de^ 
C2a(.. = 0,fc = 0) = ^y---- 

x27r(5(sinh(^3) + sinh(^4) - sinh(^2) - sinh(^i)) 

g-/3A(cosh(03)+cosh(6l4))^]^ _ g-/3A{cosh(6»3)+cosh{04)-cosh(ei)-cosh(6»2))') 

cosh(6'i) + cosh(6'2) — cosh(^^3) — cosh(6'4) 

^ Yl fa^Laiaii^s - i'^, - «vr, 61, 62)0 
010203 04 

A/f 3 

X&a^aM - tn, e, - IT,, 9,, ^3)0. (A23) 

As the 4-particle form factors are proportional to (sinh(^3) + sinh(^^4) — sinh(^2) ~ sinh(^i)) 
(the Lorentz pre-factor for the matrix element) one might believe it is immediate that this 
expression vanishes once the Fourier transform, limfc^o / dxe^^^ , is taken and so makes no 
contribution to the susceptibility. However the need to regulate the form factor leaves this 
ambiguous. Nevertheless, after the regulation C*26 ends up making no contribution to the 
susceptibility. It will however make a contribution to the NMR relaxation rate. Hence some 
of the details needed to compute C*26 will be dealt with in the context of that computation 
(see Appendix B and Section 2.C). 
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APPENDIX B: COMPUTATION OF THE CORRELATOR FOR THE NMR 

RELAXATION RATE, 1/Ti 

In order to compute 1/Ti we must evaluate the correlator, 

C{x = 0,uj = 0) = J dte''''{M^{0,t)M^{0,0)). (Bl) 

The lowest order contribution arising from the evaluation of the thermal trace takes the 
form 

Mo^(0,t)Mo^(0,0),_, order- ^^iW = / |^fle-^^™^^(^)Ee-*^(™^^(^^)-™^^(^»+^*(^-^"^-^ 

xe^^^»(A,(^)|Mi(0,0)|A,,(ei))(^ai(^i)|Mi(0,0)|A,(e)), (B2) 

where Sa is the spin of particle a. We have assumed the field, H, is aligned along the 
3-direction. Although we perform the calculation at finite H, the form-factors themselves 
retain their H = form, a feature of the model's underlying integrability. Finite H merely 
breaks the degeneracy of the triplet state with the consequent energy shifts seen above. 
For the purposes of this computation, we are interested in the regime if ^ T ^ A. This 
permits setting e^^^°- to 1, provided we are willing to tolerate errors of 0{H/T). Performing 
then the sums, J2aai ; o^^^r the different types of excitations leaves us with 

CJt) = -2 I ^^e-^^™^^(^)e-**^(™^^^(^^)-^°^^^(^«cos(m) 

xfi'sho - tn,e,)f^J{e, - t7T,e){i + o{h/^)). (B3) 

Substituting the expression for the form-factors, /23°, from Section 3 into the above, and 
then performing the necessary Fourier transform, leaves us with 

c,{u = 0) = ^ / de ; In LS ^ + ^^^/^^ + ^^^/^)^- ^^4) 

^ •> -Uinh^^e) + ^cosh(e) 



For T < A this reduces to Ci{uj = 0) ^ ^e-^^(log(4T/if) - 7), where 7 is Euler's 
constant. This is the result found in - a logarithmic dependence on H indicative of 
ballistic transport. 

The next order in the computation, essentially computing terms of (9(e~^^^), is of the 
form 

dOi d92 dO^ de^ 



x{A,M)AaM\Moi^,0)\AaM)AaM)) 



= C2l{t) + C22{t) + C2z{t) + C2,{t) + C25(t) + C2^{t). (B5) 

Here we have introduced the same notation employed to evaluate the second order contri- 
bution to the susceptibility. The definitions of €21 are the same as those in ( [A14| ) but for 
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changing Mq to Mq and shifting energies by a Zeeman term. As in the susceptibihty com- 
putation, C2i(t) is an iU-defined term proportional to 5(0), but is cancelled off by similar 
terms coming from the partition function. Similarly, C22 and C23 are disconnected terms 
related to Ci . They give a contribution of the form 

^ J ,Jsmh\9) + ^cosh{9) 

If wc were to add similarly disconnected terms coming from matrix elements with a greater 
number of particle numbers, we would find a resummation of the form: 

C2i(a; = 0) + C22{uj = 0) + (723(0; = 0) + higher order disconnected terms 

_ 2A e-^^^°^^(^)cosh^(^) 1 

'^•1 ^sinh2(^) + f cosh(^) 1 + Se-f^^^^^^'^) " ^ ^ 

This type of resummation was discussed in Section 2. A. 2. 

The remaining terms are connected. C24,{t) is given by (we again set terms of the form 

e^^^ to 1) 

C24{t) J2 [ ^^3 ^-/3A(cosh(ei)+cosh(e2))^-itA(cosh(e3)-cosh(e2)) co3(ij-^) 

4 010203 27r 27r 27r 

x|/i^°ai,oi,02(^3 - ^1 - iTT, ^i, ^2)c/s^°a3 (^2 " ^TT, ^3) 

+ E -5:^1 (^2l)^„t1 (^13)/S,a„a3 (^1 " 9, - ZTT, ^2, ^l)c/Sl " ^TT, ^3) 
040^ a2 0.^ 

+ (^2^^3)|. (B8) 

To evaluate this expression we must again regulate the four particle form-factors appearing 
in the above by removing the singularities arising when two rapidities equal one another. 
For example, we regulate the first four particle form-factor appearing in the above via 

/a3°ai,ai ,02(^3 " ^^T, ^1 - iTT, ^i, ^2)0 = 

y V'(^32 - i7r)| COSh(ei)(n V')GSi,ai,a2(^3 - ^TT, ^1 - ^1, ^2) 

+ (sinh(^2) - Sinh(^3))5-i,((n V')<^Si,ai,a2(^3 - ^1 - ZTT, ^1 - ^2)) }, (B9) 

where H = ''/'(^3i)''/'(^3i — iii + iri)ip {9 12 — i7^)4'{9i2 — if])- Regulating the other form-factors 
similarly, we find after a long computation 

C2a{ou = 0) = ^ / ci^ici^2e-^^(^°^'^(^^)+<=°^^(^^))| f^l^(«i^l^(^2) - sinh(^3)) 
256 J I V I smh(6'3) | 
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X 



g23COth^(^23/2). 



'{l2^23COsh(^i) + 12(smh(^2) - sinh(^3))(^ + ^ " 



17A 



6»3 =cosh " 1 (cosh (6*2 ) 



^13 






6*13 


<— > — 


•>) 



(BIO) 



We perform a similar procedure on C25. As with the susceptibihty, C25 must and does 
generate an identical contribution to C24. 

The remaining term to evaluate is C26- This term made no contribution to the suscep- 
tibility but does make a contribution to the relaxation rate, l/Ti- C26 takes the form 



n U\ ^ /■ C?^2 <i6'3 '^6'4 /3A(cosh(6»i)+^ 



cosh(6'2)) 



^ g-itA(cosh(6»3)+cosh(6»4)-cosh(6»i)-cosh(6»2))git-ff(So3+Sa4-Soi-Sa2) 

X /a2,Si,a4,a3(^2 - ^TT + ^£2, 6'i - ZTT + iei, 9^, 6*3)0 

TV/'"'" 

X fa3%4,ar,a2{^3 - tTl + ^63, 6*4 - ^TT + ^€4, 6'i, 92)0- 



(Bll) 



Again we must regulate this expression by discarding terms proportional to the 1/e's. To 
exhibit such terms we deform the contours 6*3 and 6*4 via 

^3 ^ ^3 + «7r; 
6*4 ^ ^4 + in. 

In doing so, we deform through a scries of poles whose residues we thus pick up. Taking 
these into account, we end up with 

C26(i) = -^ ^ f deid^2rf^4e-^^^'=°'^(^i)+'°''^^^^))e-^^*^°^'^(^^)cos(//i) 



x< e 



01020304 

iAtcosh(e2) /a2°ai,a4,a3(^2 - + ^62, - ^TT + iei, ^4, ^3) 



V^(6'i3 - in + iei) 



63=61 



x/o3°o4,oi,o2(^i - «7r + ie3, 6*4 - ivr + ie4, 6*1, 6*2) 



,i At cosh(ei ) /o2!ki,04,03(^2 - iTT + ie2, ^1 - iTT + jei, Oa, 03 



iIj{623 - iTT + ie2) 
x/64,oi,02(^2 - iTT + iea, ^4 - «7r + ^£4, 9u ^2)| 

_ J_ ^ /■ d^id^2ci^3e-''^(^°^'^^'^)+'=°^^('^))e**^^°^^('^^)cos(i/t) 

01 a'5a.'!a4 



6*3 =^2 
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f^%ua4,a3i^2 - m + ie2, 0i - in + iei, 9^, e^ + in - ie^) 



x< e 



iAt cosh(e2 ) ■/ 02,01,04,03 



■0(^14 — in + iei) 



$4=01 



X/af°a4,ai,02(^3, Oi - in + is^, Oi, 62) 



I ^itAcosh(ei) /a2?oi,o4,o3(^2 - + ^€2, 9^ - in + isi, 04, 63 + in - iea) 



6>4=6'2 



.Ml 



-0(^24 -in + ie2) 

X/sr,°a4,ai,02(^3, O2 - in + ^64, ^1, ^2)| 

+ ^ ^ /'ci^icie2ci^3ti^4e-^^('°''^('^)+'°'^('^^) 
01020304 

^gjtA(cosh(6»i)+cosh(e2)+cosh(03)+cosh(e4)) QQg(^}Jf'j 
X/a2,°ai,04,03(^2 " ^TT, ^1 - OT, 64 + in, 63 + in) fas°a4,ai,a2{^3, O4, ^1, O2) 
= EC26i(i). (B12) 



i=l 



As we are interested in C2q{uj ~ 0), we immediately see that the last three terms, C263 to 
C265, may be neglected as they are only non-zero for frequencies, lu, in excess of 2A (provided 
H < A)'. 

Focusing then upon C261, we obtain upon performing the necessary regulation 



C26i(t) = J deide2rf^4e-^''('°''^^^^)+'°'^(^^»cos(i/i)e-^^*(^°^'^(^*)-'^°^^(^^ 

x(sinh(^4) - sinh(^2))(n0)Ck,o4,o3(^2 - ^7^, 9, - tn, 64, 9i) 

1 

X I COsh(^i) (H ij)G7},a4,aua2 (^1 " ^^T, 64 - in, 9i, 62) 
^ 2 

+(sinh(^2) - sinh(^4))9-,e((n V')G'i4,oi,o2(^i - + ic, 64 - in, e,, ^2)) 

n = 0(^21)0(^24 - ^7^)V^(^2l - ^7^)0(^l4 - ^7^)0;(^4l) 
1 

n = V'(^i4 + ^e)V'(^l2 - ^7^ + k)V'(^4i - ^7^)V'(^42 - ^7^)V'(^12) (B13) 
2 

To evaluate the above expression, we first Fourier transform which then leads us to consider 
the following expressions: 



)) 



(n0)(n0)GS,W,o4,o3G™3%,oi,o2 - {e^ - -e^) = 

1 2 

(nv')G'S.i,o4,o39-..((n0')G'S.4,oi,o2) + - -Oi) 

1 ^2 ^ 

= (n0)C%,o4,o3Ck.oi,o2(5-^.(n0) + (^^ -^^) 
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^24COth^(^24/2) 



Putting everything together then yields 



[i + o{dl,) + o{el,) + o{el,)). 



C26i(a; = 0) 



16 



sinh(6'4 



(sinh(^2) - sinh(^ 



^^i4COth^(^24/2) 



X 



1 2 



+ vr^) 



x(l + 0(^2 + 0(^12)' + 



e4=cosh-i(cosh(e2)+//) 



+ (ii- ^ -H) 



127rAe-2'^^^|^(log(^) - 7)(1 + + 0{T/^)). 



The evaluation of C262 yields an identical contribution to the relaxation time. 



44 



REFERENCES 



[1] M. Takigawa, T. Asano, Y. Ajiro, M. Mekata, and Y. Uemura, Phys. Rev. Lett. 76 
(1996) 2173. 

F. D. M. Haldane, Phys. Lett. 93A (1983) 464. 
L Affleck, Phys. Rev. B 41 (1990) 6697. 
I. Affleck and R. Weston, Phys. Rev. B 45 (1992) 4667. 

E. Sorensen and L Affleck, Phys. Rev B 49 (1994) 13235. 
J. Sagi and L Affleck, Phys. Rev. B. 53 (1996) 9188; |cond-mat/9506137| . 
S. Sachdev and K. Damle, Phys. Rev. Lett. 78 (1997) 943; |cond-mat/9610115 . 
K. Damle and S. Sachdev, Phys. Rev. B 57 (1998) 8307; |3ond-mat/9711014| . 
A. B. Zamolodchikov and AL B. Zamolodchikov, Ann. of Phys. 120 (1979) 253. 
A. B. Zamolodchikov and Al. B. Zamolodchikov, Nucl. Phys. B 379 (1992) 602. 
S. Fujimoto, J. Phys. Soc. Jpn. 68 (1999) 2810; |cond-mat/99054T5| . 
S. Sachdev and K. Damle, J. Phys. Soc. Jpn. 69 (2000) 2712; |cond-mat/99U^U31 . 

G. Delfino, G. Mussardo, and P. Simonetti, Nucl. Phys. B 473 (1996) 469. 

F. A. Smirnov, Form Factors in Completely Integrable Models of Quantum Field Theory, 
World Scientific, New Jersey (1992). 

J. Cardy and G. Mussardo, Nucl. Phys. B 410 (1993) 451. 

G. Delfino and G. Mussardo, Nucl. Phys. B 455 (1995) 724. 
G. Delfino and J. Cardy, Nucl. Phys. B 519 (1998) 551; |hep-th/9712TTT . 



J. Cardy and G. Mussardo, Nucl. Phys. B 410 (1993) 451. 

A. LeClair and G. Mussardo, Nucl. Phys. B 552 (1999) 624; |hep-th/ 9902075 . 

H. Saleur, Nucl. Phys. B 567 (2000) 602; |hep-th/99090T9| . 

G. Delfino, |cond-mat/010118q . 

G. Mussardo, |hep-th/0l032T^ . 

R. M. Konik, A. LeClair, and G. Mussardo, Int. J. of Mod. Phys. A 11 (1996) 2765. 
A. LeClair, F. Lesage, S. Sachdev, and H. Saleur, Nucl. Phys. B 482 (1996) 579; |cOT 
mat/ 9606101 . 



F. Lesage and H. Saleur, Nucl. Phys. B 490 (1997) 543; |cond-mat/9612050 . 

F. Lesage and H. Saleur, Nucl. Phys. B 493 (1997) 613; |cond-mat/9611025| . 

V. Korepin, N. Boguliubov, and A. Izergin, Quantum Inverse Scattering Method and 

Correlation Functions, Cambridge University Press, Cambridge, 1993. 

H. Castella, X. Zotos, and P. Prelovsek, |cond- mat / 94 1 1 05| ; X. Zotos and P. Prelovsek, 

|cond-mat/95IMJ7^ ; X. Zotos, F. Naef, and P. Prelovsek, |cond-mat/96TTUU7 . 

M. Liischer, Nucl. Phys. B 135 (1978) 1. The relevant conserved charge here is non-local. 

For a discussion of local conserved charges see A. M. Polyakov, Phys. Lett. B 72 (1977) 

224, and Y.Y. Goldschmidt and E. Witten, Phys. Lett. B 91 (1980) 392. 

M. Garst and A. Rosch, |cond-mat/0102109| . 

S. Fujimoto and N. Kawakami, J. Phys. A 31 (1998) 465; |cond-mat/9710313| . 

J. Balog and M. Niedermaier, Nucl. Phys. B 500 (1997) 421; |hep-th/961203g . 

J. Balog, Nucl. Phys. B 419 (1994) 480. 

M. Karowski and P. Weisz, Nucl. Phys. B 139 (1978) 455. 

R. M. Konik and A. W. W. Ludwig, |cond-mat/98Tim^ . 

A. Tsvehck and P. Wiegmann, Adv. in Phys. 32 (1983) 453. 



45 



[37] A. Tsvelick, Sov. Phys. JETP 66 (1987) 221. 

[38] P. Wiegmann, JETP Lett. 41 (1985) 95. 

[39] F. A. Smirnov, Commum. Math. Phys. 132 (1990) 415. 

[40] D. Bernard and A. LeClair, Nucl. Phys. B 399 (1993) 709. 



46 



FIGURES 



0.6 



0.4 



0.2 











• TBA: Exact 
■* — < TBA: Low T 

* FF: 2+4p 




1 



T/A 



FIG. 1. Plots of the zero- field susceptibility computed both from the TBA equations and from 
the form factor expansions. The first of these is an exact numerical solution of the TBA equations 
for the 0(3) sigma model. The second is arrived from a small temperature expansion in powers of 
e~^^ of these same equations. The final plot gives the form factor computation of the susceptibility. 
We have truncated the form factor expansion at the four particle level. 
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FIG. 2. The zero-field susceptibility of a Maxwellian gas is compared here to both the exact 
susceptibility of the 0(3) NLSM and the susceptibility of the 0(3) NLSM computed via a form 
factor expansion. 



48 



0.8 



0.6 



0.4 



0.2 











A — A Resummed FF: 2+4 p 



• TBA; Exact 




1 



T/A 



FIG. 3. The zero-field susceptibility of the 0(3) NLSM as computed using a resummed form 
factor expansion is compared both with the exact result coming from the TBA equations and the 
unresummed form factor susceptibility. 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

T/A 

FIG. 4. In this plot we present the form factor computation of the Drude weight, D, of the 
spin conductance. As with the susceptibihty, both the unresummed and resummed computation 
give roughly the same answer. 
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FIG. 5. In this log-linear plot we present the form factor computation of the NMR relaxation 
rate, 1/Ti, as a function of H for a variety of temperatures. We plot a normalized rate, the ratio 
of l/Ti{H) with l/Ti{H = A/36). 
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